Astronomy & Astrophysics manuscript no. cm 


©ESO 2011 


August 22, 2011 





o 

(N 



< 



o 

u 



Mass profiles and c — M DM relation 
in X-ray luminous galaxy clusters 

S. Ettori 1 - 2 , F. Gastaldello 3 ' 5 - 6 , A. Leccardi 3 ' 4 , S. Molendi 3 , M. Rossetti 3 , D. Buote 5 , and M. Meneghetti 1,2 



INAF, Osservatorio Astronomico di Bologna, via Ranzani 1, 1-40127 Bologna, Italy 
INFN, Sezione di Bologna, viale Berti Pichat 6/2, 1-40127 Bologna, Italy 
INAF, IASF, via Bassini 15, 1-20133 Milano, Italy 

Universita degli Studi di Milano, Dip. di Fisica, via Celoria 16, 1-20133 Milano, Italy 
Department of Physics and Astronomy, University of California, Irvine, CA 92697-4575 
Occhialini Fellow 



Received 24 June 2010 / Accepted 14 September 2010 
^ ■ ABSTRACT 



Context. Galaxy clusters represent valuable cosmological probes using tests that mainly rely on measurements of cluster masses and 
. baryon fractions. X-ray observations represent one of the main tools for uncovering these quantities. 

Aims. We aim to constrain the cosmological parameters f2 m and ag using the observed distribution of the both values of the concen- 
trations and dark mass within ifeoo and of the gas mass fraction within JJsoo- 

Methods. We applied two different techniques to recover the profiles the gas and dark mass, described according to the Navarro, 
Frenk & White (1997, ApJ, 490, 493) functional form, of a sample of 44 X-ray luminous galaxy clusters observed with XMM-Newton 
in the redshift range 0.1 — 0.3. We made use of the spatially resolved spectroscopic data and of the PSF-deconvolved surface bright- 
ness and assumed that hydrostatic equilibrium holds between the intracluster medium and the gravitational potential. We evaluated 
<— j , several systematic uncertainties that affect our reconstruction of the X-ray masses. 

O . Results. We measured the concentration C2oo> the dark mass M200 and the gas mass fraction in all the objects of our sample, pro- 

viding the largest dataset of mass parameters for galaxy clusters in the redshift range 0.1 — 0.3. We confirm that a tight correlation 
between C200 and A/200 is present and in good agreement with the predictions from numerical simulations and previous observations. 
When we consider a subsample of relaxed clusters that host a low entropy core, we measure a flatter c — M relation with a total 
scatter that is lower by 40 per cent. We conclude, however, that the slope of the c — M relation cannot be reliably determined from the 
fitting over a narrow mass range as the one considered in the present work. From the distribution of the estimates of C200 and M200, 
with associated statistical (15-25%) and systematic (5-15%) errors, we used the predicted values from semi-analytic prescriptions 
calibrated through N-body numerical runs and obtain erg f2^ eo ± - 03 — 0.45 ± 0.01 (at 2a level, statistical only) for the subsample 
of the clusters where the mass reconstruction has been obtained more robustly and erg fi^ 56±0 ' 04 = 0.39 ± 0.02 for the subsample 
■ of the 11 more relaxed LEC objects. With the further constraint from the gas mass fraction distribution in our sample, we break the 

degeneracy in the erg — fi m plane and obtain the best-fit values <r 8 m 1.0±0.2 (0.83 ±0.1 when the subsample of the more relaxed 
, objects is considered) and fi m = 0.26 ± 0.02. 

Conclusions. We demonstrate that the analysis of the distribution of the C200 — M200 — /gas values represents a mature and compet- 
itive technique in the present era of precision cosmology, even though it needs more detailed analysis of the output of larger sets of 
ON ' cosmological numerical simulations to provide definitive and robust results. 

o ■ 

Key words, galaxies: cluster: general - intergalactic medium - X-ray: galaxies - cosmology: observations - dark matter. 

1. Introduction i? 2 oo = c 2 oo x r s and the cluster's volume V = 4/37ri?2oo i s 
J\. ■ equal to M2oo/(200/9 C z ), where M200 is the cluster gravitating 
5^ ■ The distribution of the total and baryonic mass in galaxy clus- mass within Rwq] with this prescription, the structural proper- 
CG ; ters is a fundamental ingredient to validate the scenario of struc- ties of DM halos from ga i axies to galaxy clusters are dependent 
ture formation in a Cold Dark Matter (CDM) Universe. Within on me halo masSi with systems a t higher masses less concen- 
this scenario, the massive virialized objects are powerful cosmo- trated Moreover, the concentration depends upon the assembly 
logical tools able to constrain the fundamental parameters of a re dshift (e.g. Bullock et al. 2001, Wechsler et al. 2002, Zhao et 
given CDM model. The A^-body simulations of structure for- al 2 003, Li et al. 2007), which happens to be later in cosmolo- 
mation in CDM models indicate that dark matter halos aggre- gies with lower matter density, O m , and lower normalization of 
gate with a typical mass density profile characterized by only the linear power spec trum on scale of 8/1- 1 Mpc, a 8 , implying 
2 parameters, the concentration c and the scale radius r s (e.g. less conce ntrated DM halos of given mass. The concentration - 
Navarro et al. 1997, hereafter NFW). The product of these two mass re i a tion, and its evolution in redshift, is therefore a strong 
quantities fixes the radius within which the mean cluster den- pre diction obtained from CDM simulations of structure forma- 
sity is 200 times the critical value at the cluster's redshift [i.e. tion and is quite sens itive to the assumed cosmological parame- 
ters (NFW; Bullock et al. 2001 ; Eke, Navarro & Steinmetz 2001 ; 
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al. 2008) and Eke et al. 2001 have provided simple and powerful 
models that match the predictions from numerical simulations 
and allow comparison with the observational measurements. 

Recent X-ray studies (Pointecouteau, Arnaud & Pratt 2005; 
Vikhlinin et al. 2006; Voigt & Fabian 2006; Zhang et al. 2006; 
Buote et al. 2007) have shown good agreement between observa- 
tional constraints at low redshift and theoretical expectations. By 
fitting 39 systems in the mass range between early-type galaxies 
up to massive galaxy clusters, Buote et al. (2007) confirm with 
high significance that the concentration decreases with increas- 
ing mass, as predicted from CDM models, and require a erg, the 
dispersion of the mass fluctuation within spheres of comoving 
radius of 8 h^ 1 Mpc, in the range 0.76 — 1.07 (99% confidence) 
definitely in contrast to the lower constraints obtained, for in- 
stance, from the analysis of the WMAP 3 years data. Since it is 
based upon a selection of the most relaxed systems, these results 
assumed a 10% upward early formation bias in the concentra- 
tion parameter for relaxed halos. Using a sample of 34 massive, 
dynamically relaxed galaxy clusters resolved with Chandra in 
the redshift range 0.06 - 0.7, Schmidt & Allen (2007) high- 
light a possible tension between the observational constraints 
and the numerical predictions, in the sense that either the relation 
is steeper than previously expected or some redshift evolution 
has to be considered. Comerford & Natarajan (2007) compiled a 
large dataset of observed cluster concentration and masses, find- 
ing a normalization higher by at least 20 per cent than the re- 
sults from simulations. In the sample, they use also strong lens- 
ing measurements of the concentration concluding that these are 
systematically larger than the ones estimated in the X-ray band, 
and 55 per cent higher, on average, than the rest of the cluster 
population. Recently, Wojtak & Lokas (2010) analyze kinematic 
data of 41 nearby (z < 0.1) relaxed objects and find a normal- 
ization of the concentration - mass relation fully consistent with 
the amplitude of the power spectrum as estimated from WMAP1 
data and within la from the constraint obtained from WMAP5. 

In this work, we use the results of the spectral analysis pre- 
sented in Leccardi & Molendi (2008) for a sample of 44 X-ray 
luminous galaxy clusters located in the redshift range 0.1 — 0.3 
with the aim to (1) recover their total and gas mass profiles, (2) 
constraining the cosmological parameters as and fl m through 
the analysis of the measured distribution of C200, M200 and bary- 
onic mass fraction in the mass range above 10 14 M Q . We note 
that this is the statistically largest sample for which this study 
has been carried on up to now between z = 0.1 and z = 0.3. 

The outline of our work is the following. In Section 2, we de- 
scribe the dataset of XMM-Newton observations used in our anal- 
ysis to recover the gas and total mass profiles with the techniques 
presented in Section 3. In Section 4, we present a detailed dis- 
cussion of the main systematic uncertainties that affect our mea- 
surements. We investigate the C200 — -M200 relation in Section 5. 
By using our measurements of C200 and M200, we constrain the 
cosmological parameters as and f2 m , breaking the degeneracy 
between these parameters by adding the further cosmological 
constraints from our estimates of the cluster baryon fraction, as 
discussed in Section 6. We summarize our results and draw the 
conclusion of the present study in Section 7. Throughout this 
work, if not otherwise stated, we plot and tabulate values esti- 
mated by assuming a Hubble constant Hq = 70h^Q km s _1 
Mpc -1 and O m = 1 — fl\ = 0.3, and quote errors at the 68.3 
per cent {la) level of confidence. 

We list here in alphabetic order, with the adopted acronyms, 
the work to which we will refer more often in the present study: 
Bullock et al. (2001 - B01); Dolag et al. (2004 - D04); Eke, 
Navarro & Steinmetz (2001 - E01); Leccardi & Molendi (2008 



- LM08); Maccio et al. (2008 - M08); Navarro, Frenk & White 
(1997 - NFW); Neto et al. (2007 - N07). 

2. The dataset 

Leccardi & Molendi (2008) have retrieved from the XMM- 
Newton archive all observations of clusters available at the end of 
May 2007 (and performed before March 2005, when the CCD6 
of EPIC-MOS 1 was switched off) and satisfying the selection 
criteria to be hot (kT > 3.3 keV), at intermediate redshift 
(0.1 < z < 0.3), and at high galactic latitude (|b| > 20°). Upper 
and lower limits to the redshift range are determined, respec- 
tively, by the cosmological dimming effect and the size of the 
EPIC field of view (15' radius). Out of 86 observations, 23 were 
excluded because they are highly affected by soft proton flares 
(see Table 1 in LM08) and have cleaned exposure time less than 
16 ks when summing MOS1 and MOS2. Furthermore, 15 ob- 
servations were excluded because they show evidence of recent 
and strong interactions (see Table 2 in LM08). The spectral anal- 
ysis of the remaining 48 exposures, for a total of 44 clusters, is 
presented in LM08 and summarized in the next subsection. In 
Table[TJ we present the list of the clusters analyzed in the present 
work. 

2. 1. Spatially resolved spectral analysis 

We use gas temperature profiles measured by LM08. A detailed 
description of how the profiles were obtained and tested against 
systematic uncertainties can be found in their paper. Here we 
briefly review some of the most important points. Unlike most 
temperature estimates the one reported in LM08 have been se- 
cured by performing background modelling rather than back- 
ground subtraction. Great care and considerable effort has gone 
into building an accurate model of the EPIC background, both in 
terms of its instrumental and cosmic components. Unfortunately 
the impossibility of performing an adequate monitoring of the 
pn instrumental background during source observation resulted 
in the exclusion of this detector from the analysis. Therefore, 
we adopt the measurements obtained from the two MOS instru- 
ments (Ml and M2, hereafter) independently in the following 
analysis. 

The impact of small errors in the background estimates on 
temperature and normalization estimates was tested both by per- 
forming Monte-Carlo simulations (a-priori tests) and by check- 
ing how results varied for different choices of key parameters 
(a-posteriori tests). The detailed analysis allowed to track sys- 
tematic errors and provide an error budget including both statis- 
tical and systematic uncertainties. 

The two profiles have been analyzed both independently and 
after they were combined as described below. Ml and M2 are 
cross-calibrated to about 5% (Mateos et al 2009). The largest dis- 
crepancy appears to be in the high energy range (above 4.5 keV), 
leading to a general tendency where M2 returns slightly softer 
spectra than Ml . Since a similar comparison between M2 and pn 
shows that the latter returns even softer spectra, the M2 experi- 
ment may be viewed as returning spectra which are intermediate 
between Ml and pn in the 0.7 — 10 keV band. As consequence 
of that, a systematic shift between the Ml and M2 temperature 
profiles is present, meaning that an higher measurements is ob- 
tained with Ml . This shift is not very sensitive to the value of the 
temperature, but instead manifests itself as a difference between 
Ml and M2 in the shape of the radial temperature profile. Using 
as reference the value of gas temperature measured with M2, we 
estimate the median deviation in the different radial bins to be 
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Table 1. Sample of the galaxy clusters. 



Cluster 


Other name 




Core CI. 


Entropy CI. 


X-ray refs. 


RXCI0003 8+0203 


Abell2700 

/lULlli. / \J\J 


0.092 


ICC 


MEC 


Pr07 Cr08 


AbelB91 1 




0.097 


NCC 


HEC 


Sn08 


AbelB827 




0.098 


ICC 


MEC 


Sn08 


RXCI0049 4-2931 


AbellS0084 


0.108 


ICC 


MEC 


Si09 


Ahell2034 




0.113 


NCC 


HEC 


Ke03 Ra07 


RXCI1516 5-0056 


Ahel12051 


0.115 


NCC 


HEC 


Pr07 O08 


RXCJ2149. 1-3041 


Abell3814 


0.118 


CC 


LEC 


Cr08, Le08 


RXCI1516 3+0005 


Ahel12050 


0.118 


NCC 


HEC 


Pr07 Cr08 


RXCJ1141.4-1216 


Abelll348 


0.119 


CC 


LEC 


Pr07 Cr08 

x i \j i j v. iuu 


RXCI1044 5-0704 


Abell1084 


0.132 


CC 


LEC 


Pr07 Cr08 

X L\J 1 , \_>X wty 


Abelll068 


RXCJ1040 7+3956 


0.138 


CC 


LEC 


Wi04 Sn08 


RXCT2218 6-3853 


Abell3856 


0.138 


NCC 


MEC 


Pr07 O08 


RXCJ0605. 8-3518 


Abell3378 


0.141 


CC 


LEC 


Pr07, Cr08, Sn08 


RXCJ0020.7-2542 


Abell22 


0.142 


NCC 


HEC 


Pr07, Cr08 


Abelll413 


RXCJ1 155.3+2324 


0.143 


ICC 


MEC 


Vi05 Ba07 Sn08 Ca09 


RXCJ2048. 1-1750 


Abell2328 


0.147 


NCC 


HEC 


Pr07 Cr08 


RXCJ0547.6-3152 


Abell3364 


0.148 


NCC 


HEC 


Pr07 Cr08 

X L \J 1 > V.1 Vy U 


Abell2204 


RXC J1632 7+0534 


0.152 


CC 


LEC 


Mo07 Sa09 


RXCI0958 3-1 103 


Abell907 


0.153 


CC 


LEC 


Vi05 Cr08 


RXCT2234 5-3744 

1.W J L*£**)^ ~J /II 


Abell3888 


0.153 


NCC 


HEC 


Cr08 


RXCI2014 8-2430 


RXCT2014 8-24 


0.161 


CC 


LEC 


Cr08 


RXCI0645 4-5413 


Abell3404 


0.167 


ICC 


MEC 


Cr08 


Ahel12218 




0.176 


NCC 


HEC 


Go04 Ra07 


Abelll689 




0.183 


ICC 


MEC 


Pe98, An04, Ca09 


Ahel1383 




0.187 


CC 


LEC 


Vi05 Ta09 7h10 


Ahel1209 




0.206 


NCC 


MEC 


Ca09 


Ahell963 




0.206 


ICC 


MEC 


Sm05 Ra07 Ca09 


Abell773 




0.217 


NCC 


HEC 


Go04 Mo07 Ta09 


Ahell1763 




0.223 


NCC 


HEC 


Dn08 Ta09 


Abell2390 




0.228 


CC 


LEC 


Vi05 Mo07 Ca09 7h10 


Abell2667 




0.230 


CC 


LEC 


Ca09 


RXCI2129 6+0005 

i.w *i i i-i _y . \y vy vy \y ^/ 




0.235 


CC 


LEC 


Ta09 7h10 


Ahel11835 




0.253 


CC 


LEC 


Mo07 7h10 


RXC10307 0-2840 


Ahel13088 


0.253 


CC 


LEC 


Fi05 7h06 


Ahel168 




0.255 


NCC 


HEC 


ZhlO 


E1455+2232 


RXCJ1457.2+2220 


0.258 


CC 


LEC 


Sn08 


RXCJ2337.6+0016 




0.273 


NCC 


HEC 


Fi05, Zh06, ZhlO 


RXCJ0303.8-7752 




0.274 


NCC 


HEC 


Zh06 


RXCJ0532.9-3701 




0.275 


CC ? 


MEC 


Fi05, Zh06 


RXCJ0232.2-4420 




0.284 


Cool core remnant ? 


MEC 


Fi05, Zh06 


ZW3146 


RBS0864 


0.291 


CC 


LEC 


Mo07 


RXCJ0043. 4-2037 


Abell2813 


0.292 


NCC 


HEC 


Zh06 


RXCJ05 16.7-5430 


AbellS0520 


0.295 


NCC 


HEC 


Zh06 


RXCJ1131.9-1955 


Abelll300 


0.307 


NCC 


HEC 


Fi05, Zh06 



Notes. We quote the name of the object, the redshift adopted and the classification based on their X-ray properties. 
{Core CI): cool cores (CC), intermediate systems (ICC) and non-cool cores (NCC). 

{Entropy CI): as in Leccardi et al. (2010), low (LEC), medium (MEC) and high (HEC) entropy cores characterizing clusters with stronger cooling 
cores and more relaxed structure (LEC), more disturbed objects (HEC) and systems with intermediate properties (MEC). 

{X-ray refs.): Baldi et al. (2007, Ba07); Cavagnolo et al. (2009, Ca09); Croston et al. (2008, Cr08); Finoguenov et al. (2005, Fi05); Govoni et al. 
(2004, Go04); Kempner et al. (2003, Ke03); Morandi et al. (2007, Mo07); Pratt et al. (2007, Pr07); Sanderson et al. (2009, Sa09); Sivanandam 
et al. (2009, Si09); Snowden et al. (2008, Sn08); Vikhlinin et al. (2005, Vi05); Wise et al. (2004, Wi04); Zhang et al. (2006, Zh06); Zhang et al. 
(2010, ZhlO). 



4.8 % in the inner bin, 8.9 % in the following 4 bins, 10 % from 
the 5th bin upwards. The two profiles are then combined by a 
weighted mean and a further systematic error is added, as de- 
scribed in Leccardi & Molendi (2008; see Section 5.3): i.e. 2, 3 
and 5 percent increases are considered at 0.3 < r/i?i 80 < 0.36, 



0.36 < r/R 18 o < 0.45 and r/R 18 o > 0.45, respectively. For 
this purpose, i?iso as defined in Tab. 3 in Leccardi & Molendi 
(2008) is considered. An error of the same amount is propagated 
in quadrature with the statistical error. 
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Fig. 1. {Left) Surface brightness profile in the 0.7 — 1.2 keV band (black filled circles) of Abelll835 compared with the profiles of 
the background components. The open diamonds show the count rate predicted from the background spectral model in the annulus 
10-12 arcmin and rescaled for the mean vignetting correction of 0.472 at those radii: the instrumental component (NXB; green), the 
photon component (CXB + galactic foregrounds; blue) and the total background (sky + instrumental; red). The dashed lines show 
the background profiles that we have used in our analysis: the "photon" background (blue), which is constant and corresponds to the 
value in the outer annulus rescaled to the center, and the instrumental background profile (green), increasing with radius in order to 
consider the over-correction of this component. The red dashed line shows the total background that we have subtracted from our 
source plus background profile, with its associated one a statistical error (red dotted lines) obtained with a Monte Carlo simulation. 
Note that the intensity of the background components and their relative contribution vary significantly from cluster to cluster. (Right) 
Example of the PSF-corrected background-subtracted surface brightness profile as obtained after the analysis outlined in Sect. 2.2. 
This example refers to Abelll835, one of the objects with the largest smearing effect due to the combination of the telescope's PSF 
and the centrally peaked intrinsic profile. 



On the other hand, no significant effect is noticed when the 
values of the normalization of the thermal model K obtained 
from the two different instruments are compared. The combined 
profile is then the direct result of the weighted mean of the two 
estimates. 

Unlike in LM08, where the focus was on the measure of the 
Tgas profile in outer regions, here we need to recover a detailed 
description of both the T gas and surface brightness Sh profiles 
at large and small radii. A significant improvement compared 
to the treatment by LM08 has been the correction of the spec- 
tral mixing between different annuli caused by the finite PSF of 
the MOS instruments. We adopted the cross-talk modification 
of the ancillary region file (ARF) generation software (using the 
crossregionarf parameter of the argen task of SAS), treating 
the cross-talk contribution to the spectrum of a given annulus 
from a nearby annulus as an additional model component (see 
Snowden et al. 2008). This is a thermal model with parameters 
linked to the thermal spectrum fitted to the nearby annulus and 
associated to the appropriate ARF file of that region (i.e., the 
usual ARF familiar to X-ray astronomers). We found the cor- 
rection to be important in particular to the first two annuli used 
in the analysis. The annuli have been therefore fitted jointly in 
XSPEC version 12 (Arnaud 1996), which allows to associate dif- 
ferent models to different RMF and ARF files. A comparison of 
the values obtained with this modelling and the values quoted in 



Snowden et al. (2008) for the 16 clusters in common with our 
sample give results in agreement within the errors. 

2.2. Analysis of the surface brightness profile 

We extend the spectral analysis presented in LM08 with a spatial 
analysis of the combined exposure-corrected M1-M2 images. 

We extract surface brightness profiles from MOS images 
in the energy band 0.7 — 1.2 keV, in order to keep the back- 
ground as low as possible with respect to the source. For this 
reason, we avoid the intense fluorescent instrumental lines of Al 
(~ 1.5 keV) and Si (~ 1.8 keV) (LM08). To correct for the 
vignetting, we divide the images by the corresponding exposure 
maps. From the surface brightness profiles, we subtract the back- 
ground that is estimated starting from the spectral modelling of 
the background components in the external ring 10-12 arcmin 
(see LM08 for details on the adopted models). We recall here 
that in the procedure of LM08 the normalizations of the back- 
ground components are the only free parameters of the fit and 
that the galactic foreground emission, the cosmic X-ray back- 
ground and the cosmic ray induced continuum give a significant 
contribution in the 0.7 — 1.2 keV energy range. The intensities 
of the background components in the annulus 10-12 arcmin are 
given by the count rates predicted by the best fit spectral model 
in this region. In order to associate errors to these count rates, 



4 



S. Ettori et al.: Mass profiles and c — Mdm relation in X-ray luminous galaxy clusters 
Abell1835, z=0.2532 
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Fig. 2. Example of the results of the two analyses adopted for the mass reconstruction. (Top and middle panels, left) Gas density 
profile as obtained from the deprojection of the surface brightness profile compared to the one recovered from the deprojection of 
the normalizations of the thermal model in the spectral analysis; observed temperature profile with overplotted the best-fit model 
(from Method 1). (Top and middle panels, right) Data (diamonds) and models (dashed lines) of the projected gas density squared 
and temperature (from Method 2). (Bottom, left) Constraints in the r s — c plane with the prediction (in green) obtained by imposing 

the relation C200 = 4.305/(1 + z) x (M 2 oo/10 14 /iroo M o) from M08 - (Bottom, right) Gas mass fraction profile obtained 

from Method 1 (gray) and Method 2 (red). 



we perform a simulation within XSPEC: we allow the normal- 
izations of the background components to vary randomly within 
their errors, we obtain the count rates associated to this fake 
model and we iterate this procedure. The error on the level of the 
background components is the width of the distribution of the 
simulated count rates. Using these values in the outer annulus, 



we reconstruct the background profile at all radii. The "photon" 
components (CXB and galactic foreground) are affected by vi- 
gnetting in the same way as the source photons and, therefore, di- 
viding by the exposure map effectively corrects also these back- 
ground components for the vignetting. In order to reconstruct 
the "photon" background profile, it is thus sufficient to rescale 
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the count rate for the mean vignetting in the outer annulus (con- 
stant blue profile in Fig. [TJ. On the contrary, the instrumental 
background does not suffer from vignetting and, therefore, di- 
viding the image by the exposure map "mis-corrects" this com- 
ponent. In order to consider this effect, we divide the correspond- 
ing count rate by the vignetting profile (that we derive from the 
exposure map in the 0.7—1.2 keV), obtaining the growing green 
curve in Fig.Q] The total background profile (red line in Fig. |T} 
is the sum of the photon (blue) and instrumental (green) profiles. 

The surface brightness profiles Sb(r) have been first ex- 
tracted from the combined images and binned by requiring a 
fixed number of 200 counts in each radial bin to preserve all the 
spatial information available. After the background subtraction, 
they have been corrected for the PSF smearing. For this purpose, 
a sum of a cusped /3— model and of a (3— model (Cavaliere & 
Fusco-Femiano 1978) with seven free parameters, f m (r) = 



v 0.5 — 303+02/2 /. o\0.5 — 3a« 

a X X 1 ~' X ^1 + X\) + 04 (1 + X2) 

(with x\ = r/ai and X2 = r ja^), is convolved with the 
predicted PSF (Ghizzardi 2001) and fitted to the observed 
profile background-subtracted Si(r) to obtain the best-fit 
convolved model f c (di;r). Finally, to correct S&(r) for the 
PSF-convolution, we apply a correction at each radius f where 
Sb{r) is measured equal to the ratio f m (fli\f)/ f c (ai]f). An 
example of the results of the procedure is shown in Fig.Q] These 
corrected profiles are, finally, used in the following analysis up 
to the radial limit, R sp , beyond which the ratio between the 
profile and the error on it (including the estimated uncertainty 
on the measurement of the background) is below 2. 



3. Estimates of the mass profiles 

We use the profiles of the spectroscopically determined ICM 
temperature and of the PSF-corrected surface brightness esti- 
mated, as described in the previous section, to recover the X-ray 
gas, the dark and the total mass profiles, under the assumptions 
of the spherical geometry distribution of the Intracluster Medium 
(ICM) and that the hydrostatic equilibrium holds between ICM 
and the underlying gravitational potential. We apply the two fol- 
lowing different methods: 

- (Method 1) This technique is described in Ettori et al. (2002) 
and has been widely used to recover the mass profiles in re- 
cent X-ray studies of both observational (e.g. Morandi et al. 
2007; Donnarummaet al. 2009, 2010) and simulated datasets 
(e.g. Rasia et al. 2006; Meneghetti et al. 2010) against which 
it has been thoroughly tested. 

We summarize here the algorithm adopted and how it uses 
the observed measurements. Starting from the X-ray surface 
brightness profile and the radially resolved spectroscopic 
temperature measurements, this method puts constraints on 
the parameters of the functional form describing the dark 
matter Mdm, defined as the total mass minus the gas mass 
(we neglect the marginal contribution from the mass in stars 
that amounts to about 10-15 % of the gas mass in mas- 
sive systems -see, e.g., discussion in Ettori et al. 2009 and 
Andreon 2010-, and is here formally included in the A/dm 
term). In the present work, we adopt a NFW profile: 

M DM (< r) = Af t ot(< r) - M gas (< r) = 4tt r 3 s p s f(x), 



200 



Ps = Pc,z- 



3 ln(l + c) -c/(l + c) 



where x — r/r s , p c _ z = 3H z /8irG is the critical density at 

1 /2 

the cluster's redshift z, H z = H x [^a + O m (l + z) 3 ] 
is the Hubble constant at redshift z for an assumed flat 
Universe (Q m + Qa = 1), and the relation -R200 = c 200 x r & 
holds. 

The two parameters (r s , C200 ) are constrained by minimizing 
a x 2 statistic defined as 



X% 



E 



data.i 



t model 



if 



(2) 



f(x) = ln(l + x) - - 



(1) 



where the sum is done over the annuli of the spectral anal- 
ysis; Tdata are the either deprojected or observed tempera- 
ture measurements obtained in the spectral analysis; T mo< x e \ 
are either the three-dimensional or projected values of the 
estimates of T gas recovered from the inversion of the hy- 
drostatic equilibrium equation (see below) for a given gas 
density and total mass profiles; er is the error on the spec- 
tral measurements. The gas density profile, n gas , is estimated 
from the geometrical deprojection (Fabian et al. 1981, Kriss 
et al. 1983, McLaughlin 1999, Buote 2000, Ettori et al. 2002) 
of either the measured X-ray surface brightness or the esti- 
mated normalization of the thermal model fitted in the spec- 
tral analysis (see Fig. |2j. In the present study, we consider 
the observed spectral values of the temperature and evaluate 
2modei by projecting the estimates of T gas over the annuli 
adopted in the spectral analysis accordingly to the recipe in 
Mazzotta et al. (2004) and using the gas density profile ob- 
tained from the deprojection of the PSF-deconvolved surface 
brightness profile (see Sect. 2.2). We exclude the deprojected 
data of the gas density within a cutoff radius of 50 kpc be- 
cause the influence of the central galaxy is expected to be not 
negligible, in particular for strong low-entropy core systems. 
The values of T gas are then obtained from 

n gas Af to t(< r ) _ d(n gas x T gas ) 
— (j(im a 5 — -, j (j) 

where G is the universal gravitational constant, m a is the 
atomic mass unit and /i=0.61 is the mean molecular weight 
in atomic mass unit. To solve this differential equation, we 
need to define a boundary condition that is here fixed to the 
value of the pressure measured in the outermost point of 
the gas density profile, P out = F gas (i? sp ) = n gas (R sp ) x 
T gas (R sp ), where T gas (i? sp ) is estimated by linear extrapo- 
lation in the logarithmic space, if required. The systematic 
uncertainties introduced by this assumption on P out are dis- 
cussed in the next section. Note that by applying Method 1 
the errors on the gas density do not propagate into the esti- 
mates of the parameters of the mass profile and are used both 
to define the range of the accepted values of P out and to eval- 
uate the uncertainties on the gas mass profiles. The allowed 
range at la of the two interesting parameters, r s and C200, is 
defined from the minimum and the maximum of the values 
that permit \t t0 be lower or equal to min(xr) + 1- The 
average error on the mass is then the mean of the upper and 
lower limit obtained at each radius from the allowed ranges 
at la- of r s and C2oo- Only for the purpose of estimating the 
profile of M gas (< r), and eventually to provide the extrapo- 
lated values, the deprojected gas density profile is fitted with 
the generic functional form described in Ettori et al. (2009) 
and adapted from the one described in Vikhlinin et al. (2006), 

% as = n g as,o [r/r Ct o) x (1 + (r/r C;0 r) x 
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Fig. 3. (F;>if 2 panels on the left) Relative errors on M200 (black) and C200 (red) estimated by the two methods. The median values 
are indicated by a dashed line. (3rd and 4th panel from left) Ratios between the best-fit result on the scale radius r s (3rd panel) and 
on -R500 (4th panel) and outermost radius reached with the spatial analysis. 



- (Method 2) The second method follows the approach de- 
scribed in Humphrey et al. (2006) and Gastaldello et al. 
(2007) where further details of this technique ( in particular 
in Appendix B of Gastaldello et al. 2007) are provided. We 
assume parametrizations for the gas density and mass pro- 
files to calculate the gas temperature assuming hydrostatic 
equilibrium, 



T( r ) = T u ngas '° - ^ maG f n gas M tot dr ^ 
°n gas (r) k B n gas (r) J ro r 2 

where n gas is the gas density, n gaSj o and To are density 
and temperature at some "reference" radius ro and ks is 
Boltzmann's constant. The n gas and T(r) profiles are fitted 
simultaneously to the data to constrain the parameters of the 
gas density and mass models. The parameters of the mass 
model are obtained from fitting the gas density and temper- 
ature data and goodness-of-fit for any mass model can be 
assessed directly from the residuals of the fit. The quality of 
the data, in particular of the temperature profile, motivated 
the use of this approach rather than the default approach 
of parametrizing the temperature and mass profiles to cal- 
culate the gas density used in Gastaldello et al. (2007). We 
projected the parametrized models of the three-dimensional 
quantities, rig as and T, and fitted these projected emission- 
weighted models to the results obtained from our analysis 
of the data projected on the sky. With respect to the paper 
cited above, the XSPEC normalization have been derived 
converting the XMM surface brightness in the 0.7 — 1.2 keV 
band using the effective area and observed projected tem- 
perature and metallicity obtained in the wider radial bins 
used for spectral extraction. The models have been integrated 
over each radial bin (rather than only evaluating at a sin- 
gle point within the bin) to provide a consistent comparison. 
We considered an NFW profile of eq. [TJ for fitting the to- 
tal mass and two models for fitting the gas density profile: 
the /3 model (Cavaliere & Fusco-Femiano 1978) a double f3 
model in which a common value of beta is assumed, and a 
cusped /3 model (Pratt & Arnaud 2002; Lewis et al. 2003). 
The last two models have been introduced to account for the 
sharply peaked surface brightness in the centers of relaxed 
X-ray systems and they provide the necessary flexibility to 
parametrize adequately the shape of the gas density profiles 
of the objects in our sample when the traditional /3 model 
fails in fitting the data. 



Hereafter, we define Ma = Mdm(< Ra) (i-e. A/200 is the 
dark matter enclosed within a sphere where the mean cluster 
overdensity in dark matter only is 200 times the critical density 
at the cluster's redshift) and / gas (< Ra) is the ratio between the 
gas mass, Af gas , and the total mass, M to t = Mum + M gas , es- 
timated within Ra, where the overdensity is here estimated by 
using the total (i.e. dark+gas) mass profile. 

The best-fit values obtained for an assumed NFW dark mat- 
ter mass profiles are quoted in Table [2] In Table [3j we present 
our estimates of R200, R500 and the gas mass fraction / gas = 
Af gas /Mt t, that is hereafter considered within i?5oo to avoid 
a problematic extrapolation of the data up to i?200- In Fig- El 
we show the relative errors provided from the two methods on 
the estimates of C200 and A/200 The distribution of the statis- 
tical uncertainties is comparable, with median values of 15- 
20% on both C200 and A/2oo with Method 1 and Method 2. 
Also the distributions of the measurements of C200 and Af2oo 
are very similar, with lst-3rd quartile range of 2.70 — 5.29 
and 4.7 - 11.1 x 10 14 A/ Q with Method 1 and 3.17 - 5.09 and 
4.3 - 9.1 x 10 14 Af Q with Method 2. 

Moreover, the two methods show a good agreement between 
the two estimates of the gas mass fraction / gas (< -R500), as 
shown in Fig. We measure a median (1st, 3rd quartile) of 
0.131(0.106,0.147), and a median relative error of 12%, with 
Method 1 and 0.124(0.108, 0.155), and a relative error of 10%, 
with Method 2. 

As shown in the last two panels of Fig. [3] we note that the 
large majority of our data is able to define a scale radius r s well 
within the radial range investigated in the spectral and spatial 
analysis, allowing a quite robust constraints of the fitted param- 
eters. 

To rely on the best estimates of the concentration and mass, 
we define in the following analysis a further subsample by col- 
lecting the clusters that satisfy the criterion that the upper value 
at la of the scale radius, as estimated from the 2 methods, is 
lower than the upper limit of the spatial extension of the detected 
X-ray emission, i.e. (r s + e r J < R sp . Imposing this condition, 
we select the 26 clusters where a more robust (i.e. with well 
defined and constrained free parameters) mass reconstruction is 
achievable. 



4. Systematics in the measurements of c 2 oo, M 2 oo 

and / gas 

The derived quantities C200, M200 and / gas (< R500) are mea- 
sured with a relative statistical error of about 20, 15 and 10 %, 
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Fig. 4. Best-fit values of r$, c, gas mass fraction / gas at i?5oo and dark mass A/dm within i?200 as obtained from Method 1 (diamonds; 
the ones including red points indicate the objects where the condition (r s + e rs ) < R sp is satisfied by both methods. R sp here plotted 
as horizontal line in the upper panel). 



respectively (see Section 3 and Figures|3]and |5j. Here, we inves- 
tigate the main uncertainties affecting our techniques that will be 
treated as systematic effects in the following analysis. 

We consider two main sources of systematic errors: (i) the 
analysis of our dataset, both for what concerns the estimates of 
the gas temperature and the reconstructed gas density profile; (ii) 
the limitations and assumptions in the techniques adopted for the 
mass reconstruction. 

In Table [4] we summarize our findings tabulated as rela- 
tive median difference with respect to the estimates obtained 
with Method 1 . Overall, we register systematic uncertainties of 



(-5,+l)% on c 20 o, (-4, +3)% on M 20 o and (-l,+4)% on 
/gas(< where these ranges represent the minimum and 

maximum estimated in the dataset investigated and quoted in 
Table E] 

4. 1. Systematics from the spectral analysis 

The ICM properties of the present dataset have been studied 
through spatially resolved spectroscopic measurements of the 
gas temperature profile and deprojected, PSF-corrected surface 
brightness profile as accessible to XMM-Newton (see Section 2). 
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Fig. 5. Estimates of M200 (left), C200 (center) and / gas (< -R500) (fight) with the two methods. (Upper panels) The color code 
indicates the objects at z < 0.15 (blue), in the range 0.15 < z < 0.25 (green) and at z > 0.25 (red). (Lower panels) Distribution of 
Low (LEC), Medium (MEC), High (HEC) Entropy Core systems. 



To assess the systematics propagated through the temper- 
ature measurements, we present the results obtained with M2 
only, i.e. before any correction introduced from the harder spec- 
tra observed with Mi (see Sect. 2.1). Overall, the systematics are 
in the order of a few per cent, with the largest offset of about 4 
per cent on the concentration and mass measurements at i?5oo 
and beyond. 

When the deprojected spectral values of the gas temperature, 
instead of the projected ones, are compared with the predictions 
from the model, we measure differences below 5% (see dataset 
labelled "T 3D "). 

On the gas density profile, we investigate the role played 
from the use of a functional form instead of the values obtained 
directly from deprojection. To this purpose, we use a revised 
form of the one introduced from Vikhlinin et al. (2006) to fit 
the gas density profile and, then, we adopt it as representative of 
the gas density profile to be put in hydrostatic equilibrium with 
the gravitational potential in equation [3] The measurements ob- 
tained are labelled "fit n gas " and show discrepancies in the order 
of 1 per cent or less. 

4.2. Systematics from the mass reconstruction methods 

With the intention to assess the the bias affecting the recon- 
structed mass values, we make use of the gas temperature and 
density profiles through two independent techniques (labelled 
Method 1 and Method 2), as described in Section 3. With respect 



to Method 1, Method 2 provides differences on A/dm that are 
lower than 10 per cent, increasing from about 1 per cent at -R200 
up to 7 per cent at -R2500 ( see Tabled. The bias on / gas remains 
stable around 3^1 per cent, suggesting that some systematics af- 
fect also the estimate of Af gas . This is due to the application of 
two different functional forms in Method 1 and Method 2 over 
a radial range that extends beyond the observational limit (see, 
e.g., Fig. 0. 

The mass reconstruction of Method 1 depends upon the 
boundary condition on the gas pressure profile. In particular, to 
solve the differential equation [3] an outer value on the pressure 
is fixed to the product of the observed estimate of the gas density 
profile at the outermost radius and an extrapolated measurement 
of the gas temperature. Using a grid of values for the pressure 
obtained from the best-fit results of the gas density and tempera- 
ture profiles, we evaluate a systematic bias on the mass of about 
3 per cent, on the gas mass fraction of 1 per cent, and on C200 of 
about 1 per cent (see dataset labelled "P ut")- 

5. The 

C200 — -^200 relstion 

In this section, we investigate the C200 — Af2oo relation. We note 
that our sample has not been selected to be representative of 
the cluster population in the given redshift range and, in the 
mean time, does not include only relaxed systems. Therefore, 
the results here presented on the C200 — M200 relation have to be 
just considered for a qualitative comparison with the predictions 
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M 200 (M Q ) M 200 (M Q ) 

Fig. 6. Data in the plane (C2001 M200) use d to constrain the cosmological parameters (ft m , <r$). The dotted lines show the predicted 
relations from Eke et al. (2001) for a given ACDM cosmological model at z — (from top to bottom: erg ~ 0.9 and erg = 0.7). 
The shaded regions show the predictions in the redshift range 0.1 — 0.3 for an assumed cosmological model in agreement with 
WMAP-1, 5 and 3 years (from the top to the bottom, respectively) from Bullock et al. (2001; after Maccio et al. 2008). The dashed 
lines indicate the best-fit range at la obtained for relaxed halos in a WMAP-5 years cosmology from Duffy et al. (2008; thin lines: 
z = 0.1, thick lines: z — 0.3). Color codes and symbols as in Fig. [5] 



from numerical simulations and to assess differences or simili- 
tude with previous work on this topic. 

As we show in Fig. fusing the measurements obtained with 
Method 1, the median relation between concentration and total 
masses for CDM halos as function of redshift is represented well 
from the analytic algorithms, as in N97, E01 and B01. These 
models relate the halo properties to the physical mechanism of 
halo formation. Considering the weak dependence of the halo 
concentrations on the mass and redshift, Dolag et al. (2004) in- 
troduced a two-parameter functional form, c = cqM b / (1 + z). 
We consider this relation in its logarithmic form and fit linearly 
to our data the expression: 

log 10 (C200 x (1 + z)) = A + B x log 10 (jT^So ) ' (5) 

A minimum in the \ 2 distribution is looked for by taking into 
account the errors on both the coordinates (we use the routine 
FITEXY in IDL). The errors are assumed to be Gaussian in 
the logarithmic space, although they are properly measured as 
Gaussian in the linear space. 

We also express our results in term of the concentration C15 
expected for a dark matter halo of 10 15 H^qMq and equal to 10 A 
once the parameters in equation [5] are used. We convert to C15 
even the results from literature obtained, for instance, at different 
overdensity, as described in the Appendix. 

We measure A w 0.6 and B systematically lower than —0.1, 
with the best-fit results obtained through Method 1 that prefer, 
with respect to Method 2, a relation with slightly higher normal- 
ization (by ^10 per cent) and flatter (by 10 — 30 per cent) distri- 
bution in mass. In both cases, a total scatter of ciog 1Q c ~ 0.13 is 
measured both in the whole sample of 44 objects, where the sta- 



tistical scatter related to the observed uncertainties is still domi- 
nant, and in the subsample of 26 selected clusters. 

When a slope B = —0.1 is assumed, as measured in numer- 
ical simulations over one order of magnitude in mass almost in- 
dependently from the underlying cosmological model (see e.g. 
Dolag et al. 2004, Maccio et al. 2008), the measured normal- 
izations of the C200 — A/200 relation fall into the range of the 
estimated values for samples of simulated clusters (see Table|5j- 

All the values of normalization and slope are confirmed, 
within the estimated errors, with both the BCES bisector method 
(as described in Akritas & Bershady 1996 and implemented 
in the routines made available from M.A. Bershady) and a 
Bayesian method that accounts for measurement errors in linear 
regression, as implemented in the IDL routine LINMIXJSRR by 
B.C. Kelly (see Kelly 2007). As we quote in Table|5J with these 
linear regression methods (and after 10 6 bootstrap resampling of 
the data in BCES), we measure a typical error that is larger by 
a factor 2 — 3 in normalization and up to 6 in the slope than the 
corresponding values obtained through the covariance matrix of 
the FITEXY method. 

These values compare well with the measurements obtained 
from numerical simulations of DM-only galaxy clusters, al- 
though these simulations sample, on average, mass ranges lower 
than the ones investigated here. Recent work from Shaw et al. 
(2006) and Maccio et al. (2008) summarize the findings. The 
slope of the relation, as previously obtained from B01 and D04, 
lies in the range (—0.160, —0.083), with a preferred value of 
about —0.1. The normalizations for low-density Universe with 
a relatively higher as, as from WMAP-1, are more in agree- 
ment with the observed constraints on, e.g., C15. For instance, 
M08 find ci 5 = 4.18, 3.41, 3.56 for relaxed objects in a back- 
ground cosmology that matches WMAP-1, 3 and 5 year data, 
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respectively^ Shaw et al. (2006) measure C15 — 4.64 using a flat 
Universe with O m = 0.3 and ct 8 = 0.95. D04 for a ACDM with 
erg = 0.9 require C15 = 4.29. All these values show the sensitiv- 
ity of the normalization to the assumed cosmology, that is further 
discussed in the section where constraints on the cosmological 
parameters (fi m — erg) will be obtained through the measured 
c — M relation. Neto et al. (2007) study the statistics of the halo 
concentrations at z = in the Millennium Simulation (with an 
underlying cosmology of f2 m = 1 — fl\ = 0.25, erg = 0.9) and 
find that a power-law with B = —0.10 and C15 = 4.33 fits fairly 
well the relation for relaxed objects, with an intrinsic logarithmic 
scatter for the most massive objects of 0.092 (see their Fig. 7). 

We note, however, that, while the normalizations we mea- 
sure for a fixed slope B = —0.1 are well in agreement with 
the results from numerical simulations, a systematic lower value 
of the slope is measured, when it is left to vary. To test the ro- 
bustness of this evidence, we have implemented Monte-Carlo 
runs using the best-fit central values estimated in N-body simu- 
lations (see Appendix B for details). With almost no dependence 
upon the input values from numerical simulations and using the 
FITEXY technique that provides the results with the most sig- 
nificant deviations from B ss —0.1, we measure in the 3 sam- 
ples here considered (i.e. all 44 objects, the selected 26 objects, 
and the only 11 LEC objects) a probability of about 0.5 (1), 20 
(42) and 26 (46) per cent, respectively, to obtain a slope lower 
than the measured 1(3)<t upper limit. These result confirm that 
the systematic uncertainties present in the measurements of the 
concentration and dark mass within -R200 are still affecting the 
sample of 44 objects, whereas they are significantly reduced in 
the selected subsamples. 

Our best-fit results are in good agreement also with previ- 
ous constraints obtained from X-ray measurements in the same 
cosmology. Pointecouteau et al. (2005) measure C15 w 4.5 and 
B = —0.04 ± 0.03 in a sample of ten nearby (z < 0.15) and 
relaxed objects observed with XMM-Newton in the temperature 
range 2 — 9 keV. Zhang et al. (2006) measure a steeper slope of 
— 1 .5±0.2, probably affected from few outliers, in the REFLEX - 
DXL sample of 13 X-ray luminous and distant (z ~ 0.3) clus- 
ters observed with XMM-Newton, that, they claim, are how- 
ever not well reproduced from a NFW profile. Voigt & Fabian 
(2006) show a good agreement with B01 results and B w —0.2 
for their estimates of 12 mass profiles of X-ray luminous ob- 
jects observed with Chandra in the redshift range 0.02 — 0.45. 
A good match with the results in D04, and within the scatter 
found in simulations, is obtained with 13 low-redshift relaxed 
systems with T gas in the range 0.7 — 9 keV as measured with 
Chandra in Vikhlinin et al. (2006). Schmidt & Allen (2007), us- 
ing Chandra observations of 34 massive relaxed galaxy clusters, 
measure B = —0.45 ± 0.12 (95% c.l.), significantly steeper 
than the value predicted from CDM simulations. Leaving free 
the redshift dependence that they estimate to be consistent with 
the (1 + z) -1 expected evolution, they measure a normaliza- 
tion C15 w 5.4 ± 0.6 (95% c.l.), definitely higher than our 
best-fit parameter. Buote et al. (2007) fit the c — M relation 
from 39 systems in the mass range 0.06 — 20 x 10 14 Af Q se- 
lected from Chandra and XMM-Newton archives to be relaxed. 
Analysing the tabulated values of the 20 galaxy clusters with 
-M200 > 10 14 M Q , that include the most massive systems from 
the XMM-Newton study of Pointecouteau et al. (2005) and the 
Chandra analysis in Vikhlinin et al. (2006), we measure B = 
-0.08 ± 0.05 and c 15 * 5.16 ± 0.36. 



We refer to Appendix A for a detailed discussion of the conversions 
adopted 



Overall, we conclude however that the slope of the c — M re- 
lation cannot be reliably determined from the fitting over a nar- 
row mass range as the one considered in the present work and 
that, once the slope is fixed to the expected value of B = —0.1, 
the normalization, with estimates of C15 in the range 3.8 — 4.6, 
agrees with results of previous observations and simulations for 
a calculations in a low density Universe. 

5. 1. The subsample of Low-Entropy-Core objects 

Following Leccardi et al. (2010), we have employed the pseudo- 
entropy ratio (er = (Tin/Tout) x (-EMiN/i2MouT)~ V 3 , where 
IN and OUT define regions within «0.05 R%so an d encircled in 
the annulus with bounding radii 0.05-0.20 Riso, respectively, 
and T and EM are the cluster temperature and emission mea- 
sure) to classify our sample of 44 galaxy clusters accordingly to 
their core properties. We identify 17 High-Entropy-Core (HEC), 
11 Medium-Entropy-Core (MEC) and 16 Low-Entropy-Core 
(LEC; see Table [Qi systems. While the MEC and HEC objects 
are progressively more disturbed (about 85 per cent of the merg- 
ing clusters are HEC) and with a core that presents less evidence 
in the literature of a temperature decrement and a peaked surface 
brightness profile (intermediate, ICC, and no cool core, NCC, 
systems), the LEC objects represent the prototype of a relaxed 
cluster with a well defined cool core (CC in Table Q]) at low en- 
tropy (see also Cavagnolo et al. 2009). These systems are pre- 
dicted from numerical simulations to have higher concentrations 
for given mass, by about 10 per cent, and lower scatter, by about 
15-20 per cent, in the c - M relation (e.g. M08, Duffy et al. 
2008). 

Out of 16, eleven LEC objects are selected under the condi- 
tion that their scale radius is within the radial coverage of our 
data. We measure their c — M relation to have slightly lower 
normalization (A w 0.5 — 0.6, C15 « 3.2 — 3.7) and flatter dis- 
tribution (B = —0.4 ±0.2) than the one observed in the selected 
subsample of 26 objects, with a dispersion around the logarith- 
mic value of the concentration of 0.08, that is about 40 per cent 
lower than the similar value observed in the latter sample. This 
is consistent in a scenario where disturbed systems have an es- 
timated concentration through the hydrostatic equilibrium equa- 
tion that is biased higher (and with larger scatter) than in relaxed 
objects up to a factor of 2 due to the action of the ICM motions 
(mainly the rotational term in the inner regions and the random 
gas term above -R500X as discussed in Lau et al. (2009; see also 
Fang et al. 2009, Meneghetti et al. 2010) for galaxy clusters ex- 
tracted from high-resolution Eulerian cosmological simulations. 

6. Cosmological constraints from the 
measurements of c2oo> M200, and / gas 

N— body simulations have provided theoretical fitting functions 
that are able to reproduce the distribution of the concentration 
parameter of the NFW density profile as function of halo mass 
and redshift (e.g. NFW, E01, B01, N07). Basically, all these 
semi-empirical prescriptions provide the expected values of the 
concentration parameter for a given set of cosmological parame- 
ters (essentially, the cosmic matter density, O m , and the normal- 
ization of the power spectrum on clusters scale, erg) for a given 
mass (the estimated cluster dark mass, A/2oo> in our case) at the 
measured redshift of the analyzed object. They assume that the 
concentration reflects the background density of the Universe 
at the formation time of a given halo. The cosmological model 
influences the concentration and virial mass because of the cos- 
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mic background density and the evolution of structure formation. 
For instance, the NFW model uses two free parameters, (/, C), 
to define the collapse redshift at which half of the final mass 
M is contained in progenitors of mass > fM, with C repre- 
senting the ratio between the characteristic overdensity and the 
mean density of the Universe at the collapse redshift. We use 
(/,C) = (0.1,3000). 

B01 assume, instead, an alternative model to improve the 
agreement between the predicted redshift dependence of the con- 
centrations and the results of the numerical simulations by using 
two free parameters, F and K, where F is still a fixed fraction 
(0.01 in our study) of a halo mass at given redshift and K in- 
dicates the concentration of the halo at the collapse redshift. K 
has to be calibrated with numerical simulations and is fixed here 
to be equal to 4 (see also Buote et al. 2007 for a detailed dis- 
cussion on the role played from the parameters F and K on the 
prediction of the concentrations as function of the background 
cosmology and halo masses). M08 have revised this model by 
assuming that the characteristic density of the halo, that in B01 
scales as (1 + z) 3 , is independent of redshift. This correction 
propagates into the growth factor of the concentration parame- 
ter that becomes shallower with respect to the mass dependence 
at masses higher than 10 13 ft- _1 M Q , permitting larger concentra- 
tions at the high-mass end than the original B01 formulation. 

The prescription in E01 defines with the only parameter C a 
(equal to 28, in our analysis, as suggested in their original work) 
the collapse redshift z c through the relation D(z c )a s(M s ) = 
CgT 1 , where D(z) is the linear growth factor, a e g is the effective 
amplitude of the linear power spectrum at z = and M s is the 
total mass within the radius at which the circular velocity of an 
NFW halo reaches its maximum and that is equal to 2.17 times 
the scale radius, r s . 

As tested in high-resolution numerical simulations (see, e.g., 
N07, M08, Duffy et al. 2008), these 3 formulations provide dif- 
ferent predictions over different mass range and redshift: for 
massive systems a z < 1, as the ones under investigation in 
the present analysis, the original B01 tends to underestimate the 
concentration at fixed halo mass; its revised version after M08 
partially compensate for this difference but still shows some ten- 
sion with numerically simulated objects (see, e.g., figure 5 in 
M08); NFW overestimates the concentration, whereas E01 pro- 
vide good estimates (see, e.g. figure 2 in Duffy et al. 2008) also 
considering its simpler and more robust formulation, being de- 
pendent upon a single parameter that does not need an indepen- 
dent calibration from simulations evolved with a given back- 
ground cosmology (note, indeed, that as pointed out in M08, 
both NFW and B01 models have normalizations that, ideally, 
have to be determined empirically for each assumed cosmology 
with a dedicated numerical simulation). 

Hereafter, we consider E01 as the model of reference and use 
the other prescriptions as estimate of the systematics affecting 
our constraints. 

In particular, to constrain the cosmological parameters of in- 
terest, a% and fi m , we calculate first the concentration C2oo,y fe = 
C20o(Mi,D, m j,(7s,k) predicted from the model investigated at 
each cluster redshift for a given grid of values in mass, Mi, cos- 
mic density parameter, f2 m j, and power spectrum normalization, 
&8,k- 

Then, we proceed with the following analysis: 



1. a new mass A^oo.j and concentration C2oo,j are estimated 
from the X-ray data for given f2 m j ; 



2. we perform a linear interpolation on the theoretical predic- 
tion of C2oo,ijfc t° associate a concentration caoojfc t° the new 
mass M200.J for given il m .j and as,k\ 

3. we evaluate the merit function xt 



= Xc( fi m,j,0- 8 ,fc) = 



(C200,i — C20Qjfc) Z 



(6) 



data/* 



"200, i 



where 6200,1 is the la uncertainty related to the measured 
C200,i and a c is the scatter intrinsic to the mean predicted 
value C2oo,jk as evaluated in Neto et al. (2007; see their 
fig. 7 and relative discussion). They estimate in the mass 
bin 10 14 25 — lO 14 75 /i~ 1 M0 a logarithmic mean value of 
the concentration parameter of 0.663, with a dispersion of 
0.092, corresponding to a relative uncertainty of 0.139. We 
take into account these estimates by associating to the ex- 
pectation of C2oo,jfc a scatter equals to 10 logC200 -j fc±,Sc , where 
e c = 0.139 x logc 20 o,jfe; 

4. a minimum in the \\ distribution, \\ min> ^ s evaluated and 
the regions encompassing x\ min + (2-3, 6.17, 11.8) are es- 
timated to constrain the best-fit values and the 1, 2, 3<r in- 
tervals in the (f2 m , a%) plane shown in Fig. [7] To represent 
the observed degeneracy in the erg — plane, we quote in 
Table [6] (and show with a dashed line in Fig. [7} the best-fit 
values of the power-law fit a% 02^ — T> obtained by fitting 
this function on a grid of values estimated, at each assigned 
Q m , the best-fit result, and associated la error, of a$. 

5. A further constraint on the Q m parameter that allows us to 
break the degeneracy in the a% — Vt m plane (as highlighted 
from the banana-shape of the likelihood contours plotted in 
Fig. IT) is provided from the gas mass fraction distribution. 
We use our estimates of / gas (< -R500) = /500 from Method 
1 quoted in Table [3] We follow the procedure described in 
Ettori et al. (2009) and assume: (i) Q, h h^ = 0.0462±0.0012 
and Hq = 70.1 ± 1.3 from the best-fit results of the joint 
analysis in Komatsu et al. (2008), (ii) a depletion parameter 
at i?5oo &500 = 0.874 ± 0.023, (iii) a contribution of cold 
baryons to the total budget / cold = 0.18(±0.05)/ gas . All the 
quoted errors are at la level. Then, we look for a minimum 



in the function xj 



x) 



E 

data.z 



[^/500,i(l + /cold)/&500 — /bar.j 



(7) 



where /bar,.? = ^b/^m.j and ef t is given from the 
sum in quadrature of all the statistical errors, namely, on 
/5001 /cold, Ho, 6500 and fib- 
6. We combine the two x 2 distribution, Xtot = Xc + Xf> an d 
plot in Fig.Qthe constraints obtained from both x\ on ty an d 



Xtot 



quoting the best-fit results in Table [6] 



7. The effect of the systematic uncertainties, assumed to be nor- 
mally distributed, is also considered by propagating them in 
quadrature to the measurements of C200 and /500, as obtained 
from the analysis summarized in Table|4] The constraints ob- 
tained after this further correction are indicated with label 
"+syst" in Table [6] 

The cosmological constraints we obtain with 3 different ana- 
lytic models (E01, B01+M08, NFW) are summarized in Table 
and likelihood contours for the model of reference E01 are plot- 
ted in Fig. [7] To represent the observed degeneracy, we con- 
strain the parameters of the power-law fit a% 57^ = V. As ex- 
pected from the properties of the prescriptions, E01 provides 
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Fig. 7. Cosmological constraints in the (fi m , og) plane obtained from equations [6] and [7] by using predictions from the model by 
Eke et al. (2001). The confidence contours at 1, 2, 3cr on 2 parameters (solid contours) are displayed. The combined likelihood with 
the probability distribution provided from the cluster gas mass fraction method is shown in red. The dashed green line indicates 
the power-law fit eg = -T- The best-fit results are quoted in Table [6] A relative logarithmic scatter of 0.139 (see Sect. |6]) is 
considered in the models. Systematic uncertainties on C200 and / gas (< -R500) as quoted in Tableware also propagated. (Left) From 
the subsample of 26 clusters satisfying the condition (r s + e rB ) < i? sp ; (center) from the subsample of the LEC objects; (right) from 
all the 44 clusters. 



constraints on erg, for given ft m , that lie between the other two, 
with 7 = 0.60 ± 0.04 and T = 0.44 ± 0.02 (at 2a level; sta- 
tistical only). We break the degeneracy of the best-fit values in 
the (as, f2 m ) plane by assuming that the cluster baryon fraction 
represents the cosmic value well. We obtain that as = 1.0 ± 0.2 
and n m = 0.27 ± 0.01 (at 2a level). When the subsample of 
1 1 LEC clusters, that are expected to be more relaxed and with 
a well-formed central cooling core, is considered, we measure 
7 = 0.55 ± 0.05, T = 0.40 ± 0.02 a s = 0.85 ± 0.16 and 
n m = 0.27 ± 0.01 (at 2cr level). 

We confirm that, assumed correct the ones measured with 
E01, NFW tends to overestimate the predicted concentrations 
and, therefore, requires lower normalization as of the power 
spectrum, whereas B01+M08 compensate with larger values of 
erg the underestimate of C200 with respect to E01. 

We assess the systematics affecting our results by compar- 
ing the cosmological constraints obtained by assuming (i) differ- 
ent algorithms to relate the cosmological models to the derived 
c — M relation, (ii) biases both in the concentration parameter 
(b c = 0.9), from the evidence in numerical simulations that re- 
laxed halos have an higher concentration by about 10 per cent 
(e.g. Duffy et al. 2008), and in the dark matter (buj = 1.1) mea- 
surements from the evidence provided from hydrodynamial sim- 
ulations that the hydrostatic equilibrium might underestimate the 
true mass by 5-20 per cent (e.g. see recent work in Meneghetti et 
al. 2010). As expected, lower concentrations and higher masses 
push the best-fit values to lower normalizations of the power 
spectrum at fixed 51 m , with an offset of about 10 per cent with 
b c = 0.9 and of few per cent bm = 1.1 and M to t- 

7. Summary and Conclusions 

We present the reconstruction of the dark and gas mass from 
the XMM-Newton observations of 44 massive X-ray luminous 
galaxy clusters in the redshift range 0.1 — 0.3. We estimate a 
dark (Aftot — Af gas ) mass within i?2oo in the range (1st and 3rd 
quartile) 4 — 10 x 10 M Q , with a concentration C200 between 
2.7 and 5.3, and a gas mass fraction within i?5oo between 0.11 
and 0.16. 



By applying the equation of the hydrostatic equilibrium to 
the spatially resolved estimates of the spectral temperature and 
normalization, we recover the underlying gravitational potential 
of the dark matter halo, assumed to be well described from a 
NFW functional form, with two independent techniques. 

Our dataset is able to resolve the temperature profiles up to 
about 0.6 — O.8-R500 and the gas density profile, obtained from 
the geometrical deprojection of the PSF-deconvolved surface 
brightness, up to a median radius of 0.9i?soo- Beyond this radial 
end, our estimates are the results of an extrapolation obtained by 
imposing a NFW profile for the total mass and different func- 
tional forms for Af gas . 

We estimate, with a relative statistical uncertainty of 15 — 
25%, the concentration C200 and the mass M200 of the dark mat- 
ter (i.e. total— gas mass) halo. We constrain the C200 — 
relation to have a normalization C15 = C200 x (1 + z) x 

(M 2O o/lO 15 Af ) ~ B of about 2.9 - 4.2 and a slope B between 
—0.3 and —0.7 (depending on the methods used to recover the 
cluster parameters and to fit the linear correlation in the logarith- 
mic space), with a relative error of about 5% and 15%, respec- 
tively. Once the slope is fixed to the expected value of B = —0.1, 
the normalization, with estimates of C15 in the range 3.8 — 4.6, 
agrees with results of previous observations and simulations for 
calculations done assuming a low density Universe. We conclude 
thus that the slope of the C200 — relation cannot be reliably 
determined from the fitting over a narrow mass range as the one 
considered in the present work, altough the steeper values mea- 
sured are not significantly in tension with the results for simu- 
lated halos when the subsamples of the most robust estimates 
are considered (see Sect. 5 and Appendix B). We measure a to- 
tal scatter in the logarithmic space of about 0.15 at fixed mass. 
This value decreases to 0.08 when the subsample of LEC clusters 
is considered, where a slightly lower normalization and flatter 
distribution is measured. This is consistent in a scenario where 
disturbed systems have an estimated concentration through the 
hydrostatic equilibrium equation that is biased higher (and with 
larger scatter) than in relaxed objects up to a factor of 2 due to 
the action of the ICM motions (see e.g. Lau et al. 2009). 
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We put constraints on the cosmological parameters (erg, f2 m ) 
by using the measurements of C200 and A/200 and by comparing 
the estimated values with the predictions tuned from numerical 
simulations of CDM universes. In doing that, we propagate the 
statistical errors (with a relative value of about 15 — 25% at la 
level) and consider the systematic uncertainties present both in 
the simulated datasets (~ 20%) and in our measurements (~ 
10%; see Table |4|. To represent the observed degeneracy, we 
constrain the parameters of the power-law fit as = T and 
obtain 7 = 0.60 ± 0.03 and T = 0.45 ± 0.02 (at 2a level) 
when the E01 formalism is adopted. Different formalisms (like 
the ones in B01, revised after M08, and NFW) induce variations 
in the best-fit parameters in the order of 20 per cent. A further 
variation of about 10 per cent occurs if a bias of the order of 10 
per cent is considered on the estimates of C200 and M2oo- 

We break the degeneracy of the best-fit values in the 
(as, fi m ) plane by assuming that the cluster baryon fraction rep- 
resents the cosmic value well. We obtain that a% = 1.0 ± 0.2 and 
O m = 0.26 ± 0.01 (at 2cr level; statistical only). 

When the subsample of 1 1 LEC clusters, that are expected 
to be more relaxed and with a well-formed central cooling core, 
is considered, we measure 7 = 0.56 ± 0.04, T = 0.39 ± 0.02 
cr 8 = 0.83 ±0.1 and O ro = 0.26 ± 0.02 (at 2cr level). 

All these estimates agree well with similar constraints ob- 
tained for an assumed low-density Universe in Buote et al. 
(2007; 0.76 < a 8 < 1.07 at 99% confidence for a ACDM model 
with fi m = 0.3) and with the results obtained by analysing the 
mass function of rich galaxy clusters [see, e.g., Wen, Han & Liu 
(2010) that summarizes recent results obtained by this cosmo- 
logical tool], showing that the study of the distribution of the 
measurements in the c — Mdm — /gas plane provides a valid 
technique already mature and competitive in the present era of 
precision cosmology. 

However, we highlight the net dependence of our results on 
the models adopted to relate the properties of a DM halo to the 
background cosmology. In this context, we urge the A— body 
community to generate cosmological simulations over a large 
box to properly predict the expected concentration associated to 
the massive (> 10 14 A/ Q ) DM halos as function of as, and 
redshift. The detailed analysis of the outputs of these datasets 
will provide the needed calibration to make this technique more 
reliable and robust. 
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Appendix A: Conversion between different 
overdensity and c - M relations 

The total mass within a given overdensity A is defined in the 
present work as 

4 



M A = -ttR 3 a Ap c 



(A.l) 



where p c z — / (8ttG) is the critical density of the Universe 
at the cluster's redshift z, Ra = ca?"s is the radius within which 
the mean cluster overdensity is A times p c z and the relation with 
the concentration ca and the scale radius r s holds by definition 
of the NFW mass profile. We assume A = 200. Hereafter, we 
refer to A as any other assumed overdensity. In case it is referred 
to the background density of the Universe, p?, = fi m . z p c . 2 , it is 



Oa(l 



zf/Hl to 



straightforward to correct A by 51 m . : 
recover the definition in equation lA.il 

To convert the tabulated values to our definition of the c — M 
relation, 



C200 



Cl5 



1 + z \ M 



M 2 , 



(A.2) 



where A/15 = 10 15 /i 70 1 Af Q , we proceed as follows: 



1. by definition, Ma/(^A) is constant and Ra/ca is fixed 
from the measurement of the scale radius. Therefore, we can 
write 

M A A/200 A 



-200 



200 



(A.3) 



2. ca and C200 are related through the assumed NFW mass den- 
sity profile 



/ C2oo\ ln(l + c A ) - c A /(l + c A ) 
V ca J ln(l + C200) - c 20 o/ (1 + C200) 



A 

200' 



(A.4) 



This function is monotonic and easily to resolve numerically 
to estimate C = ca/c2oo> that is a quantity that depends 
mostly on A and only marginally on the guessed C200, as 
shown in Fig. IA.ll For instance, for A = nSfl 1 ^ 4 ^, which 
estimate the virial overdensity predicted from the spherical 
collapse model in a flat Universe with a contribution from 
dark energy (Eke et al. 2001), C = 1.34 and 1.22 at z = 
and z = 0.3, respectively, for 51 m = 0.3, with deviations 
within 2% in the range C200 = 3 — 6. 

for a given relation ca = Crj(l + z)~ x (Ma/M*) 5 , we sub- 
stitute the above relations to obtain after simple algebrical 
operations: 



C200 



coC 38 - 1 ( A Mib M 2 , 



1 + z \200 M* Mi 5 



or 



ci 5 = c C 



3B-1 ( A M\ 



V200 AT* 



(A.5) 



(A.6) 
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Fig. A.l. Numerical solution to equation IA.4I for an assumed 
C200 in the range 3 — 8 (from the thinnest to the thickest line). 

Appendix B: Monte-Carlo realizations of the 
C200 - Af 20 o relation 

We have run Monte-Carlo (MC) simulations to test the robust- 
ness of the observed deviations in the C200 — A/200 relation de- 
scribed in Section 5. We have used as input values _the best-fit 
results (defined in the following analysis as C15 and B) obtained 
in the numerical simulations from Neto et al. (2007; see their 
equations 4 and 5) and Maccio et al. (2008; see Table Al and 
A2) and listed in Tab. [5] We have considered the results for both 
the complete sample and the relaxed objects only. To each clus- 
ter in our sample with measured mass A/200, i and redshift Zi, we 
assign the concentration C2oo,i defined as 

C200,* = 10 (l (B.l) 
k = log 10 C15 x (A/ 20 o,i/10 15 ) S /(l + 2:i) +/?xei ogc , 

where R is a random value extracted from a Gaussian distribu- 
tion and £io g c is the scatter in the log-Normal distribution mea- 
sured in the numerical simulations (~ 0.13 for samples includ- 
ing all the simulated objects and ~ 0.1 for the sample of the re- 
laxed ones; the actual values are quoted in N07 and in Table Al 
and A2 of M08). We assume that (1) our LEC objects follow 
the distribution obtained for relaxed simulated clusters; (2) all 
the remaining clusters follow the distribution estimated for the 
complete simulated halo sample; (3) the 3 samples considered 
in our analysis with 44, 26 and 1 1 clusters, respectively, are built 
considering whether each object is a LEC and/or has an upper 
limit at lcr on the scale radius lower than R sp , as discussed in 
Section 5. To be conservative in our approach, we fit equation[5] 
to the distribution in the C200 — A/200 with the FITEXY technique 
that is the one that provides the most significant deviations from 
the results obtained in numerical simulations. We repeat this pro- 
cess 10,000 times and obtain the plots shown in Fig. lB.ll for each 
of the case investigated. In Table IB. fl we summarize our finding 
from which we conclude that, while the best-fit values estimated 
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Fig. B.l. Distribution of the best-fit values of the normalization 
A and slope B after 10,000 MC realizations. The input values 
(in these plots from M08 / WMAP-5) are indicated with vertical 
solid (for all the simulated objects) and dashed (for the relaxed 
ones) lines. The red solid line represents the central value for the 
corresponding sample as quoted in Table [5] The red dotted lines 
show the 1(7 uncertainties. 



for the samples of 26 and 1 1 clusters are within the overall dis- 
tribution expected in numerical simulations, the sample of 44 
clusters provides results on the slope B that lie on the lower end 
of the distribution, as probable consequence of the uncertainties 
present both on the estimates of C200 and M200 f° r the 18 clus- 
ters that are indeed not selected for the further analysis and on 
the residual bias affecting the measurements of C200 under the 
hypothesis of the hydrostatic equilibrium (see Section 5.1 and, 
e.g., Lau et al. 2009). 
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Table 2. Results on the mass reconstruction. 



Method 1 Method 2 



Cluster 




R 


y 




C200 


M20O 


AT V. A V 


7* 




C200 




M200 


Y+ (N) 

AT V 1 ^/ 


Y (N) 




kpc 


kpc 


kpc 






10 14 M Q 




kpc 












10 14 M© 






RXCJ0003.8+0203 


605 


414 


143+1 


8.06li;f 


1.90 ± 0.23 


2.3(7) 


227 ± 77 


5 


67 


± 


1 


49 


2.66 ± 0.76 


2.7(6) 


7.1(11) 


Abell3911 


836 


754 


261 +1 ° 8 

— 59 


r CCQ+1-33 
°-° 3 -1.39 


3.88 ± 0.50 


10.2(9) 


517 ±211 


3 


18 


± 





79 


5.56 ± 1.73 


10.8(8) 


9.1(13) 


Abell3827 


792 


767 


390+Si 


4 47+0-67 


6.61 ± 0.73 


3.4(9) 


345 ± 66 


4 


82 


± 





94 


5.74 ± 0.86 


2.5(8) 


19.8(25) 


RXCJ0049.4-2931 


386 


371 


71+30 

— 19 


12.77l|;?° 


0.94 ±0.16 


5.7(6) 


104 ± 17 


9 


54 


± 





84 


1.26 ± 0.16 


6.2(5) 


7.6(6) 


Abell2034 


690 


866 




2 


4fi +0.81 


17.64 ±2.17 


7.7(9) 


436 ± 110 


4 


16 


± 





68 


7.62 ± 1.47 


9.4(8) 


62.1(22) 


RXCJ15 16.5-0056 


411 


502 


563 + ? 14 


2 


7K+0.'49 
' 0_ 06 


4.73 ± 0.42 


8.1(7) 


245 ± 139 


4 


91 


± 


2 


00 


2.22 ±0.95 


4.4(6) 


5.6(4) 


RXCJ2149. 1-3041 


663 


513 


251+S 


4 


79±°'ll 


2.21 ± 0.21 


24.6(7) 


335 ± 22 


4 


08 


± 





26 


3.28 ± 0.24 


8.0(6) 


17.2(11) 


RXCJ15 16.3+0005 


624 


514 


185+% 


7 


061H1 


2.84 ±0.41 


2.8(7) 


240 ± 92 


5 


81 


± 


1 


05 


3.49 ±0.98 


2.7(6) 


15.5(7) 


RXCJ1 141.4-1216 


676 


519 


— oO 


3 


1 c+0.19 
lu — 24 


4.88 ±0.37 


13.0(7) 


422 ± 19 


3 


50 


± 





13 


4.12 ±0.22 


3.3(6) 


6.6(6) 


RXCJ1044.5-0704 


847 


566 


286+07 

— 27 


4 


cc+0.34 
25 


2.86 ± 0.18 


13.3(7) 


388 ± 63 


3 


71 


± 





42 


3.88 ±0.58 


8.0(6) 


14.8(9) 


Abelll068 


998 


1026 


564+^1 


3 


09+0.20 
u — 22 


6.40 ± 0.48 


8.0(9) 


432 ± 8 


3 


59 


± 





05 


4.86 ±0.13 


5.0(8) 


17.4(18) 


RXCJ2218.6-3853 


764 


587 


oi " -166 


3 


16± o: fl 


8.76 ± 1.62 


9.9(7) 


524 ± 175 


3 


41 


± 





78 


7.46 ± 2.29 


7.6(6) 


1.1(7) 


RXCJ0605.8-3518 


893 


598 


369+39 


4 


in+0.34 
ru_Q 34 


4.51 ± 0.36 


17.8(7) 


380 ± 43 


4 


04 


± 





32 


4.73 ± 0.46 


9.4(6) 


16.6(8) 


RXCJ0020.7-2542 


695 


603 


473+??;; 

— 154 


4 


1 7 +l.'41 
— 1 07 


10.03 ±2.67 


17.7(7) 


599 ± 205 


3 


54 


± 





82 


12.50 ±3.79 


5.2(6) 


18.2(10) 


Abelll413 


1360 


793 


287+ 2 ,? 

' —62 


5 


OQ+0:57 
o^_0 35 


6.12 ± 0.32 


6.3(8) 


280 ± 40 


5 


84 


± 





60 


5.73 ± 0.65 


5.3(7) 


6.7(14) 


RXCJ2048. 1-1750 


806 


619 


742+?S n 

— d7U 


2 


9Q+1.63 
■^o_q 21 


5.96 ± 1.12 


2.7(7) 


1365 ± 578 


1 


37 


± 


1 


40 


8.59 ±3.01 


3.0(6) 


22.6(12) 


RXCJ0547.6-3152 


847 


624 


443+lr 3 

— 71 


4 


iot; : ™ 


7.89 ± 1.51 


4.7(7) 


445 ± 158 


3 


95 


± 





95 


7.18 ±2.20 


2.9(6) 


10.9(16) 


Abell2204 


858 


837 


816+J 37 


2 


01 +0.02 

— 28 


15.93 ± 1.19 


58.3(8) 


696 ± 30 


3 


09 


± 





09 


13.19 ±0.64 


33.7(7) 


4.4(19) 


RXCJ0958.3-1103 


1088 


639 


872+?«° 

■ ' — loo 


2 


oq+0.42 

oy_Q 39 


11.94 ± 2.02 


3.8(7) 


942 ± 223 


2 


21 


± 





37 


12.03 ± 3.28 


1.9(6) 


6.4(7) 


RXCJ2234.5-3744 


745 


640 


506+?;" 

— 22V 


4 


00+2.31 

^"-1 16 


13.42 ±4.15 


1.4(7) 


648 ± 271 


3 


52 


± 


1 


23 


15.72 ±6.11 


1.7(6) 


6.7(13) 


RXCJ2014.8-2430 


999 


878 




3 




7.56 ± 0.53 


28.6(8) 


589 ± 70 


3 


28 


± 





44 


9.57 ± 1.02 


13.1(7) 


12.0(10) 


RXCJ0645.4-5413 


1287 


904 


380lgg 5 
243+?S 

— 79 


4 


58l : ^ 


7.08 ± 1.12 


7.7(8) 


323 ±174 


5 


01 


± 


1 


76 


5.69 ± 2.14 


12.4(7) 


4.3(13) 


Abell2218 


1024 


716 


6 


9^+2.46 
u — 1 48 


4.76 ± 0.74 


11.3(7) 


404 ± 129 


4 


14 


± 





62 


6.35 ± 1.57 


10.0(6) 


6.3(10) 


Abelll689 


999 


974 


211+?q 

— 19 


8 


01 +0."64 
—0 63 


7.36 ± 0.44 


16.3(8) 


238 ± 28 


7 


51 


± 





66 


7.88 ± 0.75 


9.2(7) 


25.3(21) 


Abell383 


740 


589 


435+n 5 

— U 


3 


^vJ_0 42 


4.43 ± 0.37 


27.0(6) 


505 ± 81 


3 


18 


± 





68 


5.71 ±0.91 


11.5(5) 


4.1(4) 


Abell209 


1317 


1069 


604+? 72 , 

— loo 


3 


no+0.67 
uo_q 77 


8.60 ± 1.23 


8.8(8) 


504 ± 311 


3 


35 


± 





92 


6.77 ±2.87 


8.5(7) 


12.3(16) 


Abell963 


995 


813 


— oo 


4 


oc+0:94 
°°-0.76 


6.17 ±0.83 


5.7(7) 


233 ± 58 


6 


19 


± 


1 


06 


4.20 ±0.72 


4.6(6) 


7.7(11) 


Abell773 


977 


846 


DUO o-n 

— Zoo 


3 


97 +1.49 
z '-1.05 


10.94 ±3.12 


6.6(7) 


489 ± 188 


3 


70 


± 





61 


8.42 ±2.88 


5.1(6) 


12.0(9) 


Abelll763 


1553 


1136 


192+ 1 ;! 4 

— 49 


7 


tn+2.30 
0U -3.41 


4.25 ± 0.74 


8.1(8) 


165 ±91 


8 


23 


± 


1 


76 


3.62 ±0.77 


3.2(7) 


9.2(12) 


Abell2390 


1322 


1156 


1258+L 

— 95 


2 


0fi +0 .12 
uo -0.04 


24.71 ± 1.16 


11.6(8) 


2973 ± 13 


1 


13 


± 





03 


53.68 ± 3.48 


5.3(7) 


15.6(12) 


Abell2667 


966 


885 


993 + 4 8 


2 


94+O.O8 
z ^-0.02 


15.88 ± 0.45 


4.8(7) 


867 ± 116 


2 


43 


± 





27 


13.51 ± 1.92 


2.1(6) 


1.7(6) 


RXCJ2129.6+0005 


883 


702 


41811? 


3 


71 +0.27 
' -0.38 


5.40 ± 0.44 


3.4(6) 


350 ± 55 


4 


14 


± 





46 


4.39 ±0.61 


2.0(5) 


5.1(7) 


Abelll835 


1281 


950 


866+??, 

— 143 


2 


g 4 +0.34 
u -0.09 


17.53 ± 1.41 


10.9(7) 


1122 ± 157 


2 


14 


± 





16 


20.29 ±2.89 


9.8(6) 


6.2(12) 


RXCJ0307.0-2840 


691 


951 


611±??S 


3 


1 c+0.88 
id -0.78 


10.44 ± 2.39 


5.0(7) 


301 ± 69 


5 


18 


± 





66 


5.57 ±0.94 


5.2(6) 


2.5(5) 


Abell68 


634 


746 


834+°, 57 


2 


ec+0.82 
OO_ .06 


15.96 ±1.97 


4.4(6) 


1262 ± 234 


1 


94 


± 





25 


21.78 ±4.20 


5.0(5) 


12.6(7) 


E1455+2232 


946 


752 


214+^ 


6 


09+O.53 
oz -0.51 


3.66 ± 0.29 


2.7(6) 


210 ±35 


6 


25 


± 





54 


3.35 ±0.54 


1.3(5) 


3.9(7) 


RXCJ2337.6+0016 


803 


1004 


332t?£ 


4 


qq+3.52 

yy -2.18 


6.81 ± 1.91 


1.1(7) 


499 ± 490 


3 


52 


± 


1 


11 


8.13 ±5.31 


2.3(6) 


1.5(4) 


RXCJ0303.8-7752 


906 


1007 




1 


OC+1-04 
OJ -0.09 


13.21 ±2.33 


5.7(7) 


563 ± 525 


3 


07 


± 


1 


58 


7.82 ±5.26 


4.9(6) 


2.5(5) 


RXCJ0532.9-3701 


781 


787 


2781^° 


5 


Q7+2.43 
i " -1.82 


6.88 ± 1.83 


3.3(6) 


325 ± 195 


5 


22 


± 


1 


36 


7.39 ±4.17 


3.1(5) 


10.2(6) 


RXCJ0232.2-4420 


1099 


1032 


H721V 


1 


on+0-66 
8U -0.04 


14.28 ±1.90 


12.4(7) 


515 ± 257 


3 


15 


± 





64 


6.53 ±2.71 


13.0(6) 


3.0(6) 


ZW3146 


1210 


820 


510+H 


3 


O7+0.15 
°'-0.26 


7.79 ± 0.49 


27.5(6) 


719 ± 86 


2 


64 


± 





26 


10.48 ± 1.30 


18.1(5) 


4.3(15) 


RXCJ0043. 4-2037 


940 


823 


186+Jf 
785l 4 ^ 


7 


8O+351 


4.70 ± 1.24 


10.7(6) 


142 ± 85 


9 


16 


± 


2 


67 


3.44 ± 0.88 


7.7(5) 


2.4(7) 


RXCJ05 16.7-5430 


821 


1061 


2 


41 +2!82 
41 -0.75 


10.44 ± 2.88 


1.8(7) 


462 ± 261 


3 


49 


± 





84 


6.46 ± 1.45 


2.1(6) 


2.3(7) 


RXCJ1 131.9-1955 


1285 


1091 


707+494 

' 31 -309 


2 


43+1.16 


11.31 ±2.50 


6.9(7) 


1839 ± 629 


1 


27 


± 





84 


19.82 ±5.26 


3.0(6) 


4.1(5) 



Notes. We quote the name of the object, the upper limit of the radial range investigated in the spatial (R sp ) and spectral analysis (R xsp ), the 
best-fit values of the scale radius, the concentration parameters, M200 and minimum x with the corresponding degrees of freedom. In the case of 
Method 2, we quote two minimum \ 2 , corresponding to the minima obtained from the simultaneous fits of the temperature (xt) an d gas density 
(Xn) profiles. 
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Table 3. Estimates of R200, -R500 an d the gas mass fraction. 



Method 1 Method 2 



Cluster 








/gas 










/gas 






kpc 


kpc 






-R500 


kpc 


kpc 






R> 


iOO 


RXCJ0003. 8+0203 


1231 ± 65 


824 ± 38 





.117 




n 

U 




1360 zb 122 


899 zb 50 





.097 


4 


n 

u 


ni n 


Abell3911 


1589 ± 88 


1044 ± 41 





.146 


4- 

ZZ 


n 
u 


01 7 


1773 zb 155 


1130 ± 75 





.126 


4 


n 

u 


ni ^ 


Abell3827 


1894 ± 84 


1228 ± 45 





.140 


4 


n 
u 


012 


1823 ± 87 


1184 ± 40 





.147 


4 


n 


nns 

uuo 


RXCJ0049. 4-293 1 


980 ± 59 


666 ±37 

\J \J V — 1 — ' ) \ 


o 


.143 


4- 

ZZ 


n 

u 


u.su 


1071 zb 39 

lull — i— * ' <J 


721 zb 29 





.123 


4 


n 

u 


nns 

uuo 


Abell2034 


2491 ± 140 


1569 ± 75 





.073 


4 


n 

u. 




1957 ± 108 


1267 zb 114 





.123 


4 


n 

u 


ni r 


RXCJ15 16.5-0056 


1668 ± 65 


1039 ± 38 


o 


.105 


4 


n 

U 


nnQ 

uuy 


1309 zb 159 


845 zb 98 

L ' 1 ' ' — 1 — • J \.j 





.120 


4 


n 

u 


ni ^ 

u ±0 


RXCJ2149. 1-3041 


1298 ± 52 


846 ± 30 


o 


.131 


4 


n 

u 


097 


1452 ± 36 


942 zb 26 





.101 


± 


n 

u 


nnfi 

uuu 


RXCJ15 16.3+0005 


1416 ± 95 


940 ± 54 


o 


.142 


4 




u. 


U l 4: 


1502 ± 107 


991 zb 81 





.122 


4 


n 

u 


ni <\ 


RXCJ1 141.4-1216 


1635 ± 55 


1047 ± 31 





.086 


4 


n 
u 


.012 


1551 ± 27 


1003 ± 23 





.095 


± 


n 
u 


nn^ 

uuo 


RXC J 1044. 5 -0704 


1399 ± 35 


923 ± 19 


o 


.146 


4 


n 
u 


nnQ 


1531 zb 72 


996 zb 21 





.119 


4 


n 


nn7 

uu / 


Abelll068 


1772 ± 57 


1140 ± 31 





.091 


4 


n 
u. 


nn7 


1645 zb 13 


1061 ± 7 





.105 


± 


n 


nn9 


RXCJ2218. 6-3853 


1991 ± 159 


1275 ± 84 





.099 


4- 


n 

u 


ni s 


1900 ± 167 


1222 zb 132 





.107 


4 


n 

u 


.024 


RXCJ0605.8-3518 


1613 ± 64 


1057 ± 29 


o 


.133 


4 


n 
u 


.012 


1643 ± 49 


1071 ± 25 





.125 


± 


n 

u 


nnfi 

uuu 


RXCJ0020. 7-2542 


2023 ± 228 


1329 ± 124 


o 


.062 


4 


n 

U 


ni is 


2182 zb 200 


1415 ± 82 





.060 


4 


n 

u 


nnQ 

uuy 


Abelll413 


1837 ± 64 


1207 ± 21 


o 


.161 


4 


n 

u 


m n 


1809 ± 58 


1188 zb 28 





.167 


± 


n 

u 


nn7 

UU ( 


RXCJ2048. 1-1750 


1792 ± 155 


1110 ± 80 





.132 


4 


n 
u. 


.044 


2008 zb 269 


1187 ± 109 





.114 


4 


n 

u 


n9n 

u^u 


RXCJ0547. 6-3152 


1921 ± 161 


1251 ± 85 


o 


.105 


4 


n 

U 


uo i 


1882 zb 168 


1219 ± 81 





.116 


± 


n 

u 


ni ^ 


Abell2204 


2450 ± 79 


1549 ± 44 


o 


.115 


4 


n 
u 


nns 

UUO 


2319 ± 33 


1477 zb 47 





.126 


4 


n 


nn7 

UU I 


RXCJ0958 3-1103 


2183 ± 174 


1366 zb 87 





.086 


4 


n 
u. 


ni ^ 


2191 zb 174 


1363 zb 106 

l > J\ >> J — 1 — IV/ V f 





.087 


4 


n 

u 


ni 4 

,U 


RXCJ2234. 5-3744 


2237 ± 293 


1474 ± 164 


o 


.079 


4 




u 


nfi7 

uu 1 


2377 ± 294 


1542 ± 159 





.085 


4 


n 

u 


n9 c i 

u^o 


RXCJ2014. 8-2430 


1935 ± 56 


1245 ± 32 


o 


.136 


4 
zz 


n 


ni a 


2067 ± 70 


1323 zb 16 





.120 


4 


n 

u 


nn4 

.UU'i 


RXCJ0645.4-5413 


1919 ± 133 


1243 ± 65 


o 


.161 


4 


U 


n9n 
u.^u 


1811 zb 183 


1174 ± 83 





.177 


4 


n 

u 


.022 


Abell2218 


1671 ± 120 


1100 ± 53 

1 1 V / V ' — 1 — ' ) • > 


o 


.159 


4 
zn 


U 


ni q 


1820 ± 120 


1122 ± 66 





.154 


4 


n 

u 


ni fi 


Abelll689 


1892 ± 40 


1279 ± 24 


o 


.156 


4 


n 

u 


nns 


1946 ± 54 


1304 zb 21 





.151 


4 


n 

u 


uuo 


Abell383 


1577 ± 79 


1015 ± 39 





.121 


4. 
zz 


n 
u. 


n49 


1697 zb 100 


1090 ± 17 





.101 


4 


n 


nn^ 

uuo 


Abell209 


2006 ± 125 


1267 ± 57 





.146 


4 


n 
u. 


ni ^ 

u ±0 


1873 ± 197 


1196 ± 54 





.160 


4 


n 


ni ^ 

u 10 


Abell963 


1750 ± 95 

1 1 I ' V / — 1 — i / ' J 


1153 ± 50 





.137 


4 
zl 


n 

u. 


ni ^ 


1586 ± 74 


1049 zb 36 





.164 


4 


n 

u 


ni 1 

Ul 1 


Abell773 


2100 ± 257 


1350 ± 130 





.116 


4 


n 
u. 


041 


1959 zb 170 


1140 ± 92 





.156 


4 


n 


ni q 


Abelll763 


1644 ± 105 


1079 ± 52 





.212 


4. 

ZL 


n 


n9^ 


1575 zb 88 


1028 zb 39 

1 V./ — V.' — [_ ' J 1 / 





213 


4 


n 

u. 


ni 9 


Abell2390 


2735 ± 63 


1695 ± 36 





.108 


4 


n 
u. 


ni 


3484 ± 67 


2026 zb 57 





.079 


4 


n 


nn^ 

uuo 


Abell2667 


2374 ± 36 


1478 ± 22 





.114 


4 
zn 


n 
u. 


ni s 


2259 zb 103 


1417 zb 72 





.118 


4 


n 

u 


ni ^ 

UIO 


RXCJ2129. 6+0005 


1711 ± 60 


1099 ± 30 





.165 


ZL 


n 
u. 


ni 9 


1619 ± 63 


1042 ± 16 





.177 


4 


n 

u. 


nnfi 

.UUO 


Abelll835 


2433 ± 86 


1540 ± 46 





.120 


zb 


0. 


.012 


2539 zb 100 


1583 ± 34 





1 nq 


4 


0. 


.006 


RXC 10307 0-2840 


2030 + 199 


1302 + 103 





.105 


zb 


0. 


.017 


1 695 + 78 


1114 ± 59 





.147 


zb 


0. 


.017 


Abell68 


2293 ± 127 


1457 ± 71 





.079 


± 


0. 


.008 


2549 zb 165 


1489 ± 155 





.082 


zb 


0. 


.020 


E 1455+2232 


1484 ± 46 


980 ± 26 





.160 


zb 


0. 


.013 


1445 zb 59 


954 zb 14 





.163 


zb 


0. 


.006 


RXCJ2337.6+0016 


1779 ± 192 


1178 ± 96 





.148 


zb 


0. 


.027 


1894 zb 278 


1225 zb 173 





.141 


zb 


0. 


.033 


RXCJ0303. 8-7752 


2191 ± 179 


1347 ± 93 





.116 


zb 


0. 


.016 


1888 =b 301 


1203 zb 122 





.148 


zb 


0. 


.024 


RXCJ0532.9-3701 


1784 ± 179 


1186 ± 102 





.141 


zb 


0. 


.027 


1835 ± 233 


1207 =b 105 





.136 


± 


0. 


.019 


RXCJ0232.2-4420 


2230 ± 141 


1380 ± 71 





.123 


zb 


0. 


.013 


1798 zb 167 


1152 zb 137 





.178 


zb 





.032 


ZW3146 


1875 ± 49 


1206 ± 26 





.159 


zb 


0. 


.010 


2040 zb 77 


1293 zb 22 





.135 


zb 





.005 


RXC J0043. 4-2037 


1604 ± 157 


1068 ± 82 





.176 


± 


0. 


.032 


1472 zb 95 


982 zb 133 





.199 


zb 


0. 


.033 


RXCJ05 16.7-5430 


2029 ± 246 


1273 ± 114 





.127 


zb 





.022 


1767 zb 112 


1135 =b 57 





.157 


zb 





.013 


RXCJ1 131.9-1955 


2121 ± 206 


1325 ± 93 





.155 


zb 





.023 


2513 ± 271 


1475 ± 97 





.120 


zb 


0. 


.015 



Notes. These estimates refer to the mass models obtained with two different methods (see Table[2j and are evaluated at the ovedensities determined 
from the total (i.e. dark+gas) mass profiles. All the quoted errors are at la level. 
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S. Ettori et al.: Mass profiles and c — Mom relation in X-ray luminous galaxy clusters 
Table 4. Median deviations measured in the distribution of c 2 oo, Mdm and / gas . 



Dataset 


(C200 — C2()o)/C200 


(M DM - Mom) /Mom 


(/. 


;as /gas)/ /gas 


Method 2 


-0.013 


+0.008 




+0.036 


M2 




— 01 7 




+0 00Q 


T30 


—ft 048 






+0 094 


fit r> 




-1-0 01 1 

n^U.Lf-t -1 




+0 000 


P . 
r out 


—0.011 


+0.030 




—0.014 


at -R200 


(-0.048, +0.010) 


(-0.036, +0.030) 


(" 


0.014, +0.036) 


Method 2 


- 


-0.015 




+0.035 


M? 




018 




-^0 01 
to.u±u 


T3D 




— 04fS 




+0 09^ 


fit n 
11 L Aigas 




-1-0 01 9 




— 008 


out 




-1-0 098 




— 01 ^ 


at -R500 


— 


(-0.046, +0.028) 


(" 


0.013, +0.035) 


Method 2 




-0.073 




+0.032 


M2 




-0.013 




+0.008 


T3D 




-0.059 




+0.028 


fit Ugas 




+0.004 




+0.000 


-Pout 




+0.020 




-0.009 


at -R2500 




(-0.073, +0.020) 


(- 


0.009, +0.032) 



Notes. The deviations are measured with respect to the estimates obtained from the combined M1+M2 profile with the Method 1 for the whole 
sample of 44 clusters. Dataset: (Method 2) Method 2 is used for mass reconstruction; (M2) only the T(r) profile from M2 is used; (T^d) the 
deprojected spectral measurements of T(r) are used in Method 1 instead of the projected estimates of T modo i (see Sect. 3); (fit n gas ) a model fitted 
to the gas density profile is used in Method 1; (P ou t) the outer value of the pressure is not fixed. 
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S. Ettori et al.: Mass profiles and c — 
Table 5. Best-fit values of the C200 — M200 relation. 



Mdm relation in X-ray luminous galaxy clusters 



Dataset cis A B °"iog 10 c 



All objects (44 clusters) 



Method 1 - Weighted Mean 
Method 1 - FITEXY 
Method 1 - BCES 
Method 1 - UNMIX 

Method 2 - Weighted Mean 
Method 2 - FITEXY 
Method 2 - BCES 
Method 2 - LINMIX 



o «n +0 05 
o.ou_ 05 

3.62 ±0.07 

3.78 ±0.18 

3.79 ±0.21 



hO.03 
-0.03 



3.42: 

3.21 ±0.05 
3.51 ±0.16 
3.72 ±0.20 



0.556 ±0.006 
0.558 ±0.008 
0.577 ±0.021 
0.579 ± 0.025 

0.534 ±0.004 
0.507 ±0.006 
0.545 ±0.020 
0.571 ±0.024 



-0.1 
-0.451 ± 0.023 
-0.544 ± 0.071 
-0.444 ± 0.077 

-0.1 
-0.466 ± 0.015 
-0.612 ± 0.084 
-0.493 ± 0.067 



0.193/0.116 
0.135/0.116 
0.132/0.116 
0.132/0.116 

0.203/0.119 
0.146/0.119 
0.133/0.119 
0.131/0.119 



Selected objects (26 clusters) 



Method 1 


- Weighted Mean 


4 61 + 09 

^• ol -0.09 


0. 


.664 


±0 


.008 




-0.1 





Method 1 


- FITEXY 


4.06 ±0.17 





.608 


±0 


.018 


-0. 


.321 ±0.050 





Method 1 


-BCES 


3.84 ± 0.38 


0. 


.584 


±0 


.043 


-0. 


.586 ±0.116 


0. 


Method 1 


- LINMIX 


4.24 ± 0.47 


0. 


.628 


±0. 


.048 


-0. 


.370 ±0.125 


0. 



Method 2 


- Weighted Mean 


4 00+ 05 

^• uu -0.05 


0. 


.602 ±0.005 




-0.1 





Method 2 


- FITEXY 


2.87 ±0.13 





.458 ± 0.020 


-0. 


.576 ± 0.063 


0. 


Method 2 


-BCES 


3.18 ±0.57 





.502 ±0.078 


-0. 


.782 ± 0.259 





Method 2 


- LINMIX 


3.93 ± 0.45 





.594 ±0.050 


-0. 


.438 ±0.131 






only LEC objects (11 clusters) 



Method 1 


- Weighted Mean 


A 14+0.09 


0. 


.617 


±0. 


.009 




-0.1 


Method 1 


- FITEXY 


3.68 ±0.15 


0. 


.565 


±0 


.017 


-0. 


.297 ±0.051 


Method 1 


-BCES 


3.27 ±0.52 


0. 


.514 


±0. 


.070 


-0. 


.472 ± 0.229 


Method 1 


- LINMIX 


3.75 ± 0.45 





.574 


±0. 


.052 


-0. 


.279 ±0.150 



0.091/0.033 
0.081/0.033 
0.093/0.033 
0.081/0.033 



Method 2 


- Weighted Mean 


q OI +0.05 


0. 


.581 


±0. 


.005 




Method 2 


- FITEXY 


3.13 ±0.11 


0. 


.496 


±0 


.015 


-0 


Method 2 


-BCES 


3.21 ±0.75 





.506 


±0 


.102 


-0 


Method 2 


- LINMIX 


3.39 ±0.36 





.530 


±0. 


.046 


-0 



-0.1 



0.096/0.046 
0.080/0.046 
0.074/0.046 
0.071/0.046 



Simulations 



B01 




4.29 


0.632 




-0.102 


D04 




4.01 


0.603 




-0.130 


S06 - all, relaxed 




4.64,4.86 


0.667, 0.687 


-0. 


120, -0.160 


N07 - all, relaxed 




3.77,4.33 


0.576, 0.636 


-0. 


110,-0.100 


M08/WMAP-1- 


all, relaxed 


3.47,4.18 


0.540, 0.621 


-0. 


119, -0.104 


M08 / WMAP-3 - 


all, relaxed 


2.94, 3.41 


0.469, 0.533 


-0. 


088, -0.083 


M08 / WMAP-5 - 


all, relaxed 


2.98, 3.56 


0.474,0.551 


-0. 


110, -0.098 



Notes. The best-fit values refer to equation [5] and are obtained by using (i) the linear least-squares fitting with errors in both variables (FITEXY), 
(ii) the linear regression method BCES, (iii) a Bayesian linear regression method (LINMIX). In the last column, the total (otot — 2~2i (l/i ~ A — 

Bxi) 2 /N) and statistical (cr s t a t = J^f £ m/N) scatters are quoted, where yi = log 10 (c2oo(l ± z)), Xi = log 10 M200, £ yi is the statistical error 
on yi and TV is the number of objects. 
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S. Ettori et al.: Mass profiles and c — 
Table 6. Cosmological constraints on erg and fi m . 



Mdm relation in X-ray luminous galaxy clusters 



Model Ndata 7 F \c °8 fi m Xu 
adding / gas 



E01 


26 


0. 


.596 


± 





.030 





.449 


± 





.012 


33.4 


1 


.039 


+0. 
-0 


124 
l()(i 


0. 


25±° 


01 
.01 


79.3 


(LEC) E01 


11 


0. 


.558 


± 





.042 





.388 


± 





.018 


8.3 





.825 


+0 
-0 


114 
n.x:i 


0. 


26± 


1)2 

.01 


38.3 


(all) E01 


44 


0. 


.569 


± 





.026 





.408 


± 





.012 


64.0 





.850 


+0. 
-0 


0X7 

.056 


0. 


.281° 


01 
.01 


184.0 


B01+M08 


26 





.668 


± 





.040 





.547 


± 





.014 


33.4 


1 


.260 


+0. 
-0 


040 
07(i 


0. 


251° 


01 
.01 


79.1 


NFW 


26 


0. 


.718 


± 





.086 





.344 


± 





.036 


35.4 





.940 


+0. 

-0 


2 r ,2 
150 


0. 


251° 


01 
.01 


81.8 


E01 (6 C ) 


26 





.574 


± 





.032 





.418 


± 





.014 


36.8 





.939 


+0. 

-0 


108 
082 


0. 


251° 


01 
.01 


82.7 


E01 (b M ) 


26 


0. 


.591 


± 





.030 





.458 


± 





.012 


33.3 


1 


.003 


+0 
-0 


14, r , 
08!) 


0. 


261° 


02 

.01 


79.5 


E01 (b c , b M ) 


26 


0. 


.576 


± 





.032 





.423 


± 





.014 


36.8 





.936 


+0. 

-0 


102 
109 


0. 


■261°, 


(12 

.01 


82.9 


E01 (Mtot) 


26 


0. 


.588 


± 





.030 





.441 


± 





.012 


28.7 


1 


.006 


+0. 

-0 


116 
081 


0. 


251° 


01 
.01 


74.9 



Notes. These cosmological contraints are obtained from equations[6]and[7]corresponding to the confidence contours shown in Fig. [7] To represent 
the observed degeneracy, we quote the best-fit values of the power-law as fim = T. Errors at 2a (95.4%) level of confidence are indicated. 

Table B.l. Results of the 10,000 MC runs of the c — M relation fitted using the expression in equation[5] 



Model N obj mean (rms) A mean (rms) B -Bobs ± a -Pio-(Bobs) P3cr(B h s ) 



N07 


44 





.585(0.034) 


-0. 


.140(0. 


.106) 


-0 


.451 


± 





.023 


0.5 


1.5 


N07 


26 





.583(0.047) 


-0. 


.130(0. 


136) 


-0 


.321 


± 





.050 


14.5 


38.2 


N07 


11 





.570(0.051) 


-0. 


.134(0. 


153) 


-0 


.297 


± 





.051 


23.4 


47.4 


M08/WMAP-1 


44 





.610(0.041) 


-0. 


.095(0. 


135) 


-0 


.451 


± 





.023 


0.8 


1.8 


M08/WMAP-1 


26 





.591(0.057) 


-0. 


.154(0. 


.175) 


-0 


.321 


± 





.050 


24.6 


46.1 


M08/WMAP-1 


11 





.611(0.067) 


-0. 


.140(0. 


206) 


-0 


.297 


± 





.051 


29.9 


49.3 


M08 / WMAP-3 


44 





.524(0.040) 


-0. 


.073(0. 


130) 


-0 


.451 


± 





.023 


0.4 


0.9 


M08 / WMAP-3 


26 





.510(0.056) 


-0. 


.118(0. 


170) 


-0 


.321 


± 





.050 


18.0 


37.7 


M08 / WMAP-3 


11 





.526(0.065) 


-0. 


.107(0. 


200) 


-0 


.297 


± 





.051 


23.6 


42.6 


M08 / WMAP-5 


44 





.541(0.039) 


-0. 


.087(0. 


.125) 


-0 


.451 


± 





.023 


0.4 


1.1 


M08 / WMAP-5 


26 





.525(0.054) 


-0. 


.139(0. 


162) 


-0 


.321 


± 





.050 


20.5 


41.9 


M08 / WMAP-5 


11 





.544(0.062) 


-0. 


.125(0. 


.189) 


-0 


.297 


± 





.051 


26.0 


45.9 



Notes. B bs is the best-fit result quoted in Table[5] P\ a {B ^ s ) and (B bs) indicate the percentage of MC runs that provides an estimate of B 
lower than B bs + lcr and Bobs + 3a, respectively. 
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